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Backscattering  of  Scholte  Waves  and  Near-Grazing  Underwater 
Acoustic  Waves  by  a  Trench  at  a  Liquid/Solid  Interface 

Jacques  R.  Chamuel 


Abstract—  Ultrasonic  modeling  results  are  presented  on 
backscattering  of  broadband  pulsed  Scholte  waves  and  grazing 
underwater  acoustic  waves  from  a  trench  at  a  liquid/solid 
interface.  The  backscattered  waves  are  characterized  by 
interference  of  diffracted  waves  reflected  from  the  trench  bottom 
and  the  trailing  steep  side.  The  effects  of  trench  width  and 
depth  are  demonstrated. 


I.  Introduction 

Laboratory  ultrasonic  modeling  results  are  presented 
on  monostatic  backscattering  of  pulsed  Scholte  waves 
(also  known  as  Stoneley  waves  [1])  and  near-grazing  under¬ 
water  acoustic  waves  from  a  trench  at  a  liquid/solid  interface 
to  provide  physical  insight  into  seismoacoustic  wave  propaga¬ 
tion  in  an  ocean  with  range-dependent  boundary  topography 
(2]-(4].  The  basic  findings  pertain  to  the  general  problem 
of  wave  propagation  in  waveguides  with  outward  boundary 
deformations  [4]  such  as  in  shallow  water,  boreholes  with  frac¬ 
tures.  electromagnetic  waveguides  [5],  and  sound  barriers  [6). 
References  [7]-[9]  provide  other  examples  on  the  interaction 
of  Scholte  waves  with  topographic  features. 

II.  Results 

Experimental  studies  were  conducted  on  laboratory  ultra¬ 
sonic  models  of  a  trench  at  a  liquid/solid  interface  using  glass 
models  immersed  in  water  as  shown  in  Fig.  1 .  The  trench  was 
modeled  using  a  downstep  followed  by  an  upstep.  The  trench 
depth  is  denoted  by  h  and  its  width  d.  Broadband  pulsed 
ultrasonic  waves  were  generated  and  detected  using  a  2.54 
cm  diameter  1  MHz  compressional  transducer  connected  to 
a  Pulser  Receiver  Panametrics  Model  No.  5055PR.  Scholte 
waves  and  near-grazing  underwater  acoustic  waves  were  ob¬ 
tained  by  placing  the  transducer  at  the  water  glass  interface 
and  orienting  it  to  radiate  horizontally  parallel  to  the  interface 
tow'ards  the  trench.  The  described  method  of  producing  pulsed 
Scholte  waves  is  practical  for  liquid/solid  interfaces  were 
most  of  the  Scholte  wave  energy  is  in  the  liquid  and  the 
Scholte  wave  speed  is  only  about  0.2%  smaller  than  the  liquid 
compressional  wave  speed.  The  model  dimensions  and  results 
are  given  in  terms  of  acoustic  wavelength  in  the  water  at  1 
•MHz  based  on  a  compressional  sound  wave  speed  =  1466 

.Manuscript  received  .March  22.  1994;  revised  and  accepted  July  18.  1994. 
This  work  was  supported  by  the  Office  of  Naval  Research. 

The  author  is  with  Sonoquest  Advanced  Ultrasonics  Research  Wellesley 
Hills.  .MA  02181-0001  USA. 
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Fig.  1.  A  Schematic  diagram  of  the  experimental  configuration  showing 
trench  model  parameters  and  source/recciver  transducer. 

m/s.  The  glass  density,  compressional  and  shear  wave  speeds 
were  2.2  x  10^  kg/m^,  5775  m/s.  and  3468  m/s,  respectively. 
The  distance  x  between  the  transducer  and  the  leading  side 
of  the  trench  was  fixed  at  200A.  The  calculated  Scholte  wave 
speed  for  the  water/glass  model  was  Cs^h  =  1462  m/s. 

TTie  waveform  of  backscattered  waves  received  from  a 
model  with  d  =  20A,  and  /i  =  8.3A  is  displayed  in  Fig.  2. 
The  first  received  small  pulse  "A"  is  the  reflected  component 
from  the  leading  steep  side  of  the  trench.  The  second  large 
pulse  "B"  is  the  reflection  from  the  trailing  steep  side  of 
the  trench.  The  third  small  pulse  “C”  is  a  backscattered 
component  reflecting  from  both  the  trench  horizontal  bottom 
and  trailing  steep  side.  The  relative  amplitude  of  pulses  B 
and  C  depends  on  d  and  h.  Fig.  3(a)  displays  the  measured 
relative  amplitude  of  the  backscattered  waves  plotted  versus 
d/X  for  three  different  values  of  h  (3.9A.6.5A  and  8.3A).  No 
corrections  were  made  on  the  amplitude  of  the  backscattered 
waves  to  account  for  attenuation  and  geometrical  spreading. 
Each  of  the  three  curves  of  Fig.  3(a)  has  an  inflection. 
The  inflection  location  moves  to  greater  values  of  d/X  as 
the  trench  depth  is  increased.  The  inflection  is  caused  by 
destructive  interference  between  pulses  B  and  C.  The  time 
delay  between  pulses  B  and  C  was  calculated  from  the 
approximate  expression: 

T  =  2(  ^dy9  +  /i2/4  -  s/hTjm^/C^.. 

The  first  term  in  the  numerator  of  (1)  is  the  distance  from 
PI  the  top  comer  of  the  leading  side  of  the  trench  to 
“P2"  the  reflection  point  along  the  bottom  of  the  trench.  The 
second  term  is  the  distance  from  “P2”  to  “P3”  the  midpoint 
along  the  vertical  trailing  side  of  the  trench.  The  third  term 
is  the  distance  from  PI  to  P3.  Interference  between  pulses 
B  and  C  occurs  when  T  is  less  than  the  2.4  ^s  duration 
of  the  pulsed  waveform.  Fig.  4(a)  presents  the  calculated 
time  delay  between  pulses  B  and  C  plotted  versus  d/X  for 
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Fig.  2.  Waveform  of  received  backscaitered  waves  with  d  =  20A  and 
h  =  S.3A  obtained  from  contoured  air/waier/glass  model. 
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Fig.  3.  Measured  relative  amplitude  of  received  backscattcrcd  waves  from 
trench  plotted  versus  d/X  in  (a)  and  plotted  versus  dX/h-  in  (b).  The  results 
are  not  corrected  for  attenuation  or  geometrical  spreading. 

h  =  3.9A.  6.5A  and  8.3A.  TTie  intersection  of  each  curve  with 
the  horizontal  2.4  line,  representing  the  waveform  duration, 
occurs  at  rf/A  =  12.  35,  and  58  nearly  matching  the  inflection 
locations  shown  in  Fig.  3(a).  The  experimental  backscatter 
amplitude  results  of  Fig.  3(a)  are  replotted  versus  dX/h^  in  Fig. 
3(b)  showing  that  the  inflection  occurs  at  about  rfA/Zi^  =  0.9 
for  all  three  curves.  For  the  given  range  of  parameters  and  for 
T  <  5  /iS,  the  three  time  delay  curves  of  Fig.  4(a)  become 
almost  superimposed  on  each  other  when  replotied  versus 
dX/h^  as  illustrated  in  Fig.  4(b). 

The  observed  large  backscattered  waves  are  mainly  due 
to  having  a  penetrable  elastic  boundary  that  allows  seismo* 


0  50  100  150  200 

d/X 


(a) 


dX/h- 

(b) 

Fig.  4.  Calculated  time  delay  between  pulses  B  and  C  plotted  versus  d/X  in 
(a)  and  plotted  versus  dX/h^  in  (b).  The  2.4  horizontal  line  corresponds 
to  waveform  duration. 
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Fig.  5.  Waveform  of  received  backscattered  waves  with  d  =  20A  and 
h  —  8. 3 A  obtained  from  contoured  air/watcr/air  pressure>release  model. 


acoustic  interface  wave  conversion  insonifying  the  shadow 
zone.  Fig.  5  demonstrates  the  relatively  small  backscattered 
diffracted  acoustic  wave  obtained  when  the  solid  glass  is 
replaced  with  a  contoured  pressure-release  boundary.  A  con¬ 
toured  pressure-release  boundary  was  achieved  by  placing  a 
390  /im  thick  plastic  membrane  between  the  water  and  the 
glass  to  trap  a  thin  air  layer  following  the  trench  profile  to 
decouple  the  acoustic  waves  in  the  water  from  the  solid  and 
create  a  contoured  “free"  pressure-release  boundary.  The  time 
of  arrival  of  the  backscattered  diffracted  acoustic  wave  shown 
in  Fig.  5  is  approximately  equal  to  the  time  of  arrival  of  pulse 
B  (Fig.  2). 

III.  Conclusion 

Experimental  results  have  been  presented  on  backscattering 
of  Scholie  waves  and  near-grazing  underwater  acoustic  waves 
from  a  trench  at  a  liquid/solid  interface.  TTie  backscattered 
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waves  are  characterized  by  interference  of  waves  reflecting 
from  the  trench  bottom  and  the  distant  steep  side. 

The  results  indicate  that  significant  backscatter  may  be 
obtained  from  an  object  located  in  a  shadow  zone  when 
penetrable  boundaries  are  involved.  The  findings  relate  to 
interpretation  of  underwater  acoustic  wave  data  as  applied 
to  remote  mapping  of  ocean  bottom  bathymetry,  transmission 
loss  prediction,  acoustic  imaging,  and  nondestructive  testing. 
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SECTION  2 


2  EXPLANATION  FOR  ACOUSTIC  BACKSCATTERING  FROM  CANADA  BASIN 
NEAR  73.2°N,  139°W  WHERE  NO  KNOWN  BATHYMETRIC  HIGHS  EXIST 


Low-frequency  acoustic  backscattering  from  underwater  topographic  features  has  been 
used  to  characterize  ocean  boundaries  over  long  ranges.  A  survey  of  low-frequency  acoustic 
backscattering  experiments  was  reported  in  1986  by  Baggeroer  and  Dyer  [2.1].  Spindel  and 
Heirtzler  |2.2]  used  directional  towed  array  receivers  near  the  Bahamas  and  correlated  the 
acoustic  backscatter  with  topographic  features.  In  the  Arctic  Ocean,  Dyer,  Baggeroer,  Zittel, 
and  Williams  [2.3-2.51  used  two-dimensional  hydrophone  arrays  to  identify  individual 
bathymetric  features  in  azimuth  and  range,  and  superbly  correlated  low-frequency  acoustic 
backscattering  data  with  known  topographic  highs  such  as  mid-oceanic  ridges,  escarpments, 
and  continental  margins  for  ranges  over  1000  km.  Jiang  [2.6|  presented  other  low-frequency 
acoustic  backscattering  data  obtained  from  the  eastern  Arctic  Ocean. 

The  Canada  Basin  Acoustic  Reverberation  Experiment  (CANBARX)  results  reported 
by  Dyer  et  al.  [2.3]  also  generated  new  localized  acoustic  backscatter  which  were  interpreted 
as  originating  from  a  region  near  73°N  latitude,  139°W  longitude  in  an  area  of  the  Canada 
Basin  with  no  previously  known  bathymetric  highs  [2.3,  2.7,  2.8j.  Grantz  and  Hart  [2.7] 
charted  the  anomalous  acoustic  backscatter  as  a  cluster  of  three  isolated  features  (at  about 
73.2°N,  139°W,  at  72.6°N,  145. 5°W,  and  at  73.4°,  144°W)  as  shown  in  Fig.  2.1  and 
pointed  out  that  the  backscatter  is  not  produced  by  bathymetric  highs  in  the  vicinity.  In  this 
report,  the  large  central  anomalous  acoustic  feature  shown  in  Fig.  2.1  will  be  referred  to  as 
"FI",  and  the  second  largest  feature  as  "F2".  Grantz  and  Hart  obtained  new  seismic 
reflection  data  in  1993  [2.9-2.  lOj  indicating  the  absence  of  large  bathymetric  highs  or  geologic 
features  in  the  vicinity  of  73°N,  139°W  (Fig.  2.2).  Three  segments  of  the  seismic  reflection 
data  (Line  9318)  are  displayed  in  Fig.  2.2.  The  beginning  and  end  of  Line  9318  are  marked 
on  the  bottom  map. 

An  explanation  is  needed  to  comprehend  the  origin  of  the  puzzling  localized  anomalous 
backscatter  observed  by  Dyer  et  al.  [2.3|.  Chamuel  [2. 1 1|  hypothesized  that  focused  multiple 
scattering  can  explain  the  anomalous  backscatter  observed  by  Dyer  et  al.  as  explained  below. 

Two-dimensional  hydrophone  array  data  provide  azimuth  and  range  information.  Range 
is  calculated  from  the  round  trip  propagation  time  along  a  straight  line,  and  azimuth  angle  is 
determined  from  the  observed  angle  0,  of  the  returning  acoustic  waves  incident  on  the  array. 
Range  can  be  overestimated  when  acoustic  waves  propagate  first  in  different  directions  towards 
a  remote  major  topographic  feature  that  focuses  the  reflected  waves  onto  a  much  closer 
scatterer  that  in  turn  reflects  the  large  amplitude  focused  waves  back  towards  the 
source/receiver  with  the  azimuth  angle  0,.  When  focused  multiple  scattering  are  not  taken  into 
consideration,  late-arriving  backscattered  waves  approaching  the  receiver  from  the  azimuth 
direction  0,  can  be  misinterpreted  as  originating  from  a  remote  large  scatterer  at  the  same 
azimuth  angle  0,  . 
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Generalized  baihymccric  map  of  the  wcsicm  Arctic  Ocean  show- 
ing  location  of  acoustic  backscanering  features  [Dyer  et  at.  (1982)  and 
their  spatial  relation  to  (he  continental  ntargin.  Shaded  area  shows  iso* 
backscancr  levels  of  -  36  dB  and  higher.  Contours  show  iso-backscancr 
levels  of  -  36  dB  and  -  30  dB.  Isobaths,  in  kilometers,  from  Johnson  et 
aMi982|. 


Fig.  2. 1  Diagram  copied  from  A.  Grantz  and  Hart  [J.  Geophys.  Res.  89(C2), 
2105(1984)]  depicting  cluster  of  three  observed  anomalous  acoustic 
backscattering  features  near  73 °N  latitude,  139°W  longitude  where  no  known 
bathymetric  highs  exist. 
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Fig.  2.2  Seismic  reflection  profile  (line  9318)  by  Hart  and  Grantz  (1994)  demonstrating 
absence  of  bathymetric  or  geologic  features  in  the  vicinity  that  could  explain 
the  localized  cluster  of  anomalous  acoustic  backscatter  (Dyer  et  al.  1982].  The 
map  shows  the  beginning  and  end  of  line  9318  (Grantz  1994). 
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Fig.  2.3  shows  two  general  bathymetry  maps  of  the  Beaufort  Sea.  The  anomalous 
acoustic  backscattering  feature  FI  (range  353  km)  is  marked  on  the  top  map,  and  feature  F2 
(range  275  km)  is  marked  on  the  bottom  map.  According  to  Chamuel’s  hypothesis  [2. 1 1],  FI 
is  caused  by  scattering  from  the  near-elliptical  boundary  B1  (stretching  from  71.32°N, 
148. 5°W  to  7I.86°N,  152.63°W  ,  Fig.  2.3(a))  focused  on  a  bathymetric  feature  A1  (73.5°N, 
140. 4°W)  located  at  a  much  closer  range  (210  km)  than  FI  as  shown.  Similarly,  F2  is  caused 
by  focused  multiple  scattering  between  another  near-elliptical  boundary  B2  (Fig.  2.3(b))  and 
feature  A2  located  at  a  much  closer  distance  (192  km)  than  F2.  The  range  of  these  features 
can  be  overestimated  by  43-68%  when  focused  multiple  scattering  is  ignored. 


Fig.  2.3  Anomalous  acoustic  backscatter  FI  (where  no  known  bathymetric  highs  exist) 
may  be  the  result  of  scattering  by  a  major  near-elliptical  boundary  B1  focused 
onto  a  much  closer  bottom  feature  Al.  Both  S/Bl/Al/R  and  S/Fl/R  wave 
components  have  the  same  total  round  trip  travel  time  and  the  same  return 
azimuth  angle.  A  similar  hypothesis  can  be  used  to  explain  the  second 
observed  anomalous  acoustic  backscatter  F2  shown  in  the  bottom  map,  where 
the  waves  are  focused  by  B2  onto  a  second  closer  bottom  feature  A2. 
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Fig,  2.4  Ultrasonic  modeling  results  (Scale  1:1,200,000)  comparing  the  relative 
amplitude  and  time  of  arrival  of  direct  backscattered  wave  component  S/A/R 
and  late-arriving  focused  multiple  scattering  component  S/B/A/R. 


Scaled  laboratory  ultrasonic  modeling  experiments  were  conducted  to  verify  Chamuel’s 
hypothesis  and  provide  an  explanation  for  the  anomalous  acoustic  backscatter  observed  by 
Dyer  et  al.  .  The  author  chose  the  simplest  scaled  ultrasonic  model  possible  to  clearly  isolate 
the  dominant  phenomena  involved  and  compare  direct  backscattered  waves  and  multiple 
backscattered  acoustic  waves  focused  by  a  large  elliptical  reflector  segment  on  the  scatterer. 
Fig.  2.4  presents  the  experimental  setup  configuration  and  the  results  obtained  from  the  scaled 
ultrasonic  model  (horizontal  scale  1:1,200,000).  The  ultrasonic  model  consisted  of  a)  a  2-D 
water  waveguide  ,  b)  a  small  scatterer  "A"  positioned  in  the  water  at  one  focal  point  of  an 
ellipse  in  the  plane  of  the  waveguide,  c)  a  source  and  directional  receiver  "S/R"  placed  at  the 
second  focal  point  of  the  ellipse,  and  d)  a  segment  of  an  elliptical  reflector  "B"  approximately 
modeling  the  near-elliptical  bathymetry  profile  between  71.32°N,  148. 5°W  to  71.86°N, 
152.63°W  positioned  as  shown  in  Fig.  2.4.  The  elliptical  reflector  was  an  air-backed  vertical 
boundary  intersecting  the  entire  cross  section  of  the  water  waveguide.  The  receiver  was 
oriented  towards  the  scatterer  and  the  source  (3  mm  diameter)  was  excited  with  a  single 
broadband  pulse.  The  received  signal  was  amplified,  filtered,  and  displayed  on  a  digital 
oscilloscope.  Fig.  2.4. (a)  displays  the  direct  backscattered  wave  "S/A/R"  obtained  without 
the  elliptical  reflector.  The  horizontal  time  scale  was  100  us/div. ,  and  the  water  compressional 
wave  velocity  was  1480  m/s.  The  waveform  shown  in  Fig.  2.4(b)  was  acquired  with  the 
elliptical  reflector  in  place  as  shown.  The  amplitude  of  the  direct  reflection  "S/A/R"  was 
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small  compared  to  the  late-arriving  focused  backscattered  component  "S/B/A/R".  The 
propagation  time  for  S/A/R  and  S/B/A/R  was  277  us  and  457  us  respectively.  The  position 
and  orientation  of  the  elliptical  reflector  followed  the  geometry  of  Fig.  2. 3. (a)  simulating  the 
B1  boundary.  If  one  attempts  to  interpret  the  waveform  (Fig.  2.4(b))  without  taking  focused 
multiple  scattering  into  consideration,  the  large  backscattered  wave  S/B/A/R  can  be 
misinterpreted  as  originating  from  direct  backscattering  from  a  remote  large  scatterer  "F"  at 
the  same  azimuth  angle  as  diagrammed  in  Fig.  2.4. 

Another  experiment  was  conducted  using  an  irregular  near-elliptical  segment  reflector 
following  the  1500  m-depth  contour  derived  from  the  Arctic  Ocean  Bathymetry  map  [2.12] 
(Fig.  2.5).  The  amplitude  of  the  late-arriving  backscattered  wave  focused  by  the  irregular 
reflector  was  still  substantially  greater  than  the  direct  backscattered  wave.  The  results  of  this 
experiment  also  suggest  that  in  view  of  the  averaging  effect  of  the  processor  time  window  used 
by  Dyer  et  al.,  variation  in  the  position  of  source/receiver  due  to  drift  of  the  ice  camp  for 
example  would  have  a  moderate  effect  on  the  corresponding  backscattered  focused  waves.  A 
time  window  of  20s  is  equivalent  to  about  15  km  variation  in  the  reflecting  boundary  profile. 


BACKSCATTER  WAVEFORM 


WITHOUT 

ELUPTICAL 

REFLECTOR 


WITH 

ELLIPTICAL 

REFLECTOR 


1500  m 


S/A/R  S.MVA/R 


Fig.  2.5  Ultrasonic  modeling  results  similar  to  previous  figure  obtained  using  an 
irregular  near-elliptical  reflector  following  the  actual  1500  m  depth  map  profile. 
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One  important  question  to  be  answered  is:  If  Chamuel’s  hypothesis  is  correct,  then, 
scatterers  A1  and  A2  predicted  by  the  hypothesis  should  have  produced  detectable  direct 
backscattering  signals  in  the  original  Arctic  field  data.  The  author  examined  the  acoustic 
backscattering  data  in  References  12.1, 2. 3-2. 5,  2.7,  2.8)  searching  for  clues  on  the  existence 
of  the  hypothesized  ocean  bottom  features  at  A1  and  A2.  The  author  was  gratified  to  find  in 
one  of  the  9Hz  figures  in  Williams’  thesis  |2.5|  two  backscattering  signals  originating  from 
73.3°N,  142. 8°W  and  72.6°N,  145. 5°W  corresponding  precisely  to  locations  A1  and  A2 
predicted  by  Chamuel’s  hypothesis  as  shown  in  Fig.  2.6. 

Ultrasonic  modeling  results  were  presented  indicating  that  focused  multiple  scattering 
should  not  be  disregarded  in  interpreting  two-dimensional  hydrophone  array  underwater 
acoustic  data.  Future  bathymetric  measurements  need  to  be  conducted  near  73.5°N,  140. 4°W 
and  72.9°N,  147. 4°W  to  verify  the  existence  of  the  scattering  features  predicted  by  the 
proposed  hypothesis. 


Fig.  2.6  Acoustic  backscatter  data  with  lowest  contour  =  -  39  db  (copied  from  fig.  4.1 
by  Williams  1981]  revealing  the  existence  of  two  backscattering  features 
precisely  at  the  locations  A1  and  A2  as  predicted  by  Chamuel’s  focused 
multiple  scattering  hypothesis  [1994]. 
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SECTION  3 


3  EFFECT  OF  A  THIN  LOW-VELOCITY  LAYER  UNDERNEATH  A  FLOATING 
PLATE  ON  FLEXURAL  WAVE  DISPERSION 

3.1  Abstract 

Theoretical  and  experimental  results  from  a  layered  ultrasonic  model  are  presented 
providing  an  explanation  for  previously  reported  10-40%  discrepancies  between  theoretical  and 
measured  Arctic  field  results  by  Hunkins  |3. 1  j  on  the  dispersion  of  flexural  waves  in  sea  ice 
between  50-500  Hz.  The  presence  of  a  thin  low-velocity  fluid  layer  underneath  a  floating 
plate  can  substantially  decrease  the  flexural  wave  group  velocity  maximum  of  the  Scholte 
wave  branch  and  the  velocity  of  higher  frequency  components  (3.2).  Theoretical  calculations 
indicate  that,  within  a  nominal  range  of  sea  ice  parameters,  the  observed  deviation  of  Arctic 
field  data  flexural  wave  dispersion  is  due  to  prevailing  effects  of  the  compressional  velocity 
and  thickness  of  this  low  velocity  layer. 

3.2  Introduction 


Sea  ice  is  a  complex  heterogeneous  structure  |3.3-3.41.  Hunkins’  seismic  studies  on 
sea  ice  |3, 1 1  revealed  two  discrepancies  between  measured  and  calculated  flexural  wave  group 
velocity  dispersion.  The  first  discrepancy  observed  by  Hunkins  was  that  sea  ice  thickness 
determined  from  flexural  wave  dispersion  is  generally  lower  than  the  actual  thickness.  The 
second  discrepancy  in  Hunkins  field  results  is  originating  from  flexural  wave  components  to 
the  right  of  (meaning  higher  frequencies)  having  10-40  %  lower  velocities  than  predicted 
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tig.  12.  ComporiioD  oi  observed  and  thoorctlcal  diaper»ioo  for  record  J23,  April  2-f,  1058. 
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Fig.  3.1  Hunkins’  results  on  flexural  wave  dispersion  in  floating  sea  ice 

[data  copied  from  J.  Geophys.  Res.  65(10),  3469  (1960)|. 
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Fig.  3.2  Effect  of  water  loading  and  shallow  surface  cracks  on  decreasing  the  flexural  wave 
velocity  .  Model  consisted  of  Plexiglas  plate  (thickness  H  =  2.26  mm)  on  water. 
Average  crack  spacing  s  =  5  mm,  and  crack  depth  g  =  H/3. 


Previous  ultrasonic  modeling  studies  [3. 5-3. 6]  provided  an  explanation  for  the  first 
discrepancy  by  demonstrating  that  the  presence  of  a  large  number  of  shallow  cracks  in  a 
floating  plate  reduces  the  flexural  wave  group  velocity  and  shifts  the  corresponding  dispersion 
curves  towards  the  behavior  of  and  "effective"  thinner  plate.  Fig.  3.2  illustrates  the  decrease 
in  flexural  wave  phase  velocity  by  multiple  shallow  surface  cracks  and  water  loading  in  a 
floating  Plexiglas  plate. 

The  second  discrepancy  on  flexural  waves  in  Hunkins’  results  cannot  be  resolved  by 
simply  changing  the  ice  plate  thickness  to  match  the  theoretical  and  field  results.  The  author 
believes  that  by  adding  a  thin  low  velocity  layer  underneath  Hunkins’  sea  ice  plate  model  a 
much  closer  agreement  between  the  measured  and  predicted  flexural  wave  results  can  be 
obtained.  Bogorodskii  et  al.  [3.7]  reported  that  during  the  rapid  growth  of  the  ice  cover,  there 
are  two  longitudinal  wave  velocity  minima  across  the  thickness  of  a  snow-covered  ice  pack. 
One  minimum  is  at  the  snow  cover  with  a  velocity  less  than  the  velocity  of  sound  in  air,  and 
the  second  minimum  is  at  the  "clear-ice"  layer  at  the  ice/water  interface  with  a  velocity  of 
about  1050  m/s  which  is  less  than  the  sound  velocity  in  water. 

3.3  Theory 

Kheisin  [3.8]  studied  guided  wave  propagation  due  to  a  discontinuity  in  the  ocean 
sound  wave  velocity  underneath  the  ice  cover,  where  a  water  layer  about  twice  as  thick  as  the 
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ice  plate  with  a  compressional  wave  velocity  of  1450  m/s  resting  on  a  water  half-space  with 
a  velocity  of  1500  m/s.  Kheisin’s  paper  focused  on  determining  the  required  minimum  layer 
thickness  at  which  unattenuated  waveguide  propagation  is  still  possible  for  a  given  ice  plate 
thickness.  The  current  paper  is  concerned  with  the  water/ice  thin  layer  transition  zone  (about 
one  tenth  the  ice  plate  thickness)  where  the  longitudinal  wave  velocity  may  be  lower  than  the 
water  compressional  wave  velocity  [3.7].  The  shear  wave  speed  is  assumed  to  be  zero  in  the 
low-velocity  layer. 

The  theoretical  results  presented  in  this  paper  are  based  on  expanding  Kheisin’s 
characteristic  determinant  Dj  from  reference  [3.8]  for  an  isotropic  solid  plate  floating  on  a 
layered  fluid  half-space  with  a  low-velocity  upper  layer. 

Upon  expanding  Kheisin’s  determinant  D2  the  following  characteristic  equation  was  obtained: 


-  N  M  [tanht})  coth6  +  tanhS  cothc})]  +  N'^ 
M  [cothtf)  +tanh4)]  -  N  [coth6+  tanhS] 


where, 
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e  =  {TifH/V)  [l-V2/a2]°-5 


ijf  =  (2ti fd/V)  [l-V^/v^]°-^ 


f  =  frequency,  H  =  plate  thickness,  V  =  phase  velocity,  a  =  plate  compressional  wave 
velocity,  /3  =  plate  shear  wave  velocity,  p,  =  plate  density,  d  =  layer  thickness,  v  =  layer 
sound  wave  velocity,  p,  =  layer  density,  c  =  compressional  wave  velocity  of  liquid  half¬ 
space,  and  =  density  of  liquid  half-space. 

The  group  velocity  U  was  calculated  from  the  fundamental  relation 
U  =  V  +  k  dV/dk 
where  k  is  the  wavenumber. 


3.4  Results 

In  order  to  investigate  the  effect  of  a  low  velocity  layer  underneath  a  floating  plate  on 
flexural  waves,  laboratory  experiments  were  conducted  on  a  layered  ultrasonic  model  as 
described  below. 

Fig.  3.3(a)  shows  a  diagram  of  the  ultrasonic  model.  A  2.3  mm  thick  glass  plate  was 
placed  on  top  of  a  water  "half-space".  Source  and  receiver  transducers  were  placed  on  the 
free  surface  of  the  glass  plate  to  generate  and  detect  flexural  waves.  The  source  and  receiver 
transducers  were  piezoelectric  broadband  compressional  transducers  (Sonoquest  3  mm 
diameter)  placed  on  the  top  surface  of  the  glass  plate  separated  by  the  distance  x  =  15  cm. 
The  source  was  excited  with  a  single  broadband  pulse  (30  KHz-1  MHz)  and  the  received 
waveform  was  recorded  on  a  digital  oscilloscope  (Tektronix  7D20).  The  plate  dimensions  and 
water  depth  were  selected  large  enough  to  avoid  interference  from  edge  reflections.  The 
general  experimental  setup  was  previously  described  elsewhere  13.5,3.8).  The  waveform  of 
the  received  dispersive  flexural  wave  is  presented  in  Fig.  3.3(b). 

The  experiment  was  repeated  with  an  alcohol  layer  between  the  glass  plate  and  water 
(Fig.  3.3  (c)).  A  0.1  mm  thick  acetate  membrane  was  used  to  separate  the  alcohol  from  the 
water.  The  edges  of  the  membrane  were  bonded  to  the  bottom  surface  of  the  glass  plate 
forming  a  sealed  chamber  filled  with  alcohol.  The  alcohol  was  introduced  through  a  1  mm 
diameter  hole  predrilled  in  the  glass.  The  hole  was  located  away  from  the  direct  wave  path. 
The  alcohol  volume  was  controlled  to  produce  the  desired  layer  thickness.  The  alcohol  layer 
thickness  "d"  was  0.7  mm.  The  waveform  of  the  received  flexural  wave  with  the  alcohol 
layer  is  shown  in  Fig.  3.3(d).  The  alcohol  layer  slowed  down  the  flexural  wave.  The  effect 
of  the  membrane  on  the  flexural  wave  dispersion  was  negligible.  The  measured 
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Fig,  3.3  Ultrasonic  modeling  results.  Top:  waveform  of  flexural  wave  from  glass/water  model. 

Bottom;  waveform  from  glass/alcohol/water  model.  The  effect  of  the  0.1  mm  thick 
acetate  membrane  used  to  separate  the  alcohol  from  the  water  is  negligible. 


FREQUENCY  (KHz) 


Fig.  3.4  Effect  of  a  thin  low-velocity  layer  underneath  a  floating  plate  on  flexural  wave 
group  velocity  dispersion.  Theoretical  results  (solid  lines).  Experimental 
results  from  layered  model  with  d  =  0.7  mm  (open  circles). 
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Fig.  3.5  Effect  of  varying  the  thickness  d  of  the  low-velocity  layer  on  the  flexural  wave 
group  velocity.  The  model  parameters  are:  a  =  3500  m/s,  /3  =  1856  m/s, 
H  =  2  m,  p.  =  0.89  x  10'  kg/m',  v  =  1005  m/s,  p.  =  0.8  x  10'  kg/m',  c  = 
1440  m/s,  and  p,„  =  1  x  10'  kg/m'.  The  open  circles  are  from  Hunkins’  sea 
ice  flexural  wave  data. 


compressional  wave  speeds  of  the  glass  plate,  alcohol,  and  water  were  5830,  1207,  and  1480 
m/s,  respectively.  The  shear  wave  speed  of  the  plate  was  3413  m/s.  The  densities  in  the 
same  order  were  2.25,  0.8,  and  1.0  xlO'  Kg/m'. 

The  flexural  wave  group  velocity  experimental  results  are  compared  with  theoretical 
dispersion  curves  in  Fig.  3.4.  The  experimental  and  theoretical  results  agree  reasonably  well. 
The  low-frequency  components  are  practically  unaffected  by  the  presence  of  the  thin  alcohol 
layer.  The  group  velocity  of  the  higher  frequency  components  decreased  similar  to  Hunkins’ 
sea  ice  flexural  wave  field  data. 

Fig.  3.5  displays  three  theoretical  sea  ice  flexural  wave  group  velocity  curves  plotted 
as  function  of  frequency  for  three  low-velocity  layer  thicknesses  d  =  0,  0.2,  and  0.4  m.  The 
open  circles  results  in  FIG.  3  were  obtained  from  Hunkins’  Arctic  field  data  [3.1]. The  top 
curve  is  similar  to  the  theoretical  curve  provided  by  Hunkins.  One  notes  that  the  theoretical 
curve  (d=0)  deviates  by  almost  40  %  from  the  field  results  to  the  right  of  the  group  velocity 
maximum.  The  plot  dimensions  and  scales  were  selected  to  match  Hunkins’  plots.  The 
theoretical  curves  were  based  on  the  following  parameters:  a  =3500  m/s,  jS  =  1856  m/s,  p, 
=  0.89  X  10',  H  =  2  m,  v  =  1005  m/s,  =  0.8  x  10',  d  =  0.2  m,  c  =  1440  m/s,  and 

=  1  X  10'. 
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Other  theoretical  results  were  calculated  to  demonstrate  the  effect  of  varying  the 
compressional  wave  speed  V  and  density  of  the  low-velocity  layer  (Fig.  3.6  -  Fig.  3.7). 
Variations  in  V  or  d  produced  similar  effects  on  the  flexural  wave.  Nominal  changes  in 
had  little  effect  on  the  flexural  wave  dispersion  (Fig.  3.7).  For  comparison,  Fig.  3.8  displays 
the  effect  of  the  ice  shear  wave  speed  on  the  flexural  wave  group  velocity  of  the  floating  plate. 

3.5  Conclusions 

TTieoretical  and  laboratory  experimental  results  presented  in  this  paper  demonstrated 
the  effect  of  a  thin  low- velocity  layer  underneath  a  floating  plate  on  reducing  the  group 
velocity  of  low-frequency  flexural  waves  and  provided  an  explanation  for  previously  reported 
discrepancies  on  flexural  wave  dispersion  in  sea  ice  (3.1|  between  50-500  Hz.  The  results 
demonstrated  that  a  15  cm  layer  under  a  2  m  thick  ice  plate  can  reduce  the  group  velocity  of 
flexural  wave  components  by  as  much  as  40%  although  the  flexural  wavelength  is  greater  than 
the  ice  plate  thickness.  Parametric  studies  revealed  that  the  trend  observed  in  Hunkins  flexural 
wave  field  data  was  matched  by  varying  the  thickness  of  the  low-velocity  layer  and  its 
compressional  wave  velocity.  The  findings  contribute  to  fundamental  understanding  of  high 
latitude  underwater  acoustic  wave  phenomena. 
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Fig.  3.6  Effect  of  varying  the  compressional  wave  velocity  of  the  low-velocity  layer  on 
the  flexural  wave  group  velocity.  The  model  parameters  are:  a  =  3500  m/s, 
P  =  1856  m/s,  H  =  2  m,  Pj  =  0.89  x  10’  kg/m’,  d  =  0.2  m,  =  0.8  x  10’ 
kg/m’,  c  =  1440  m/s,  and  =  I  x  10’  kg/m’. 
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Fig.  3.7  Effect  of  varying  the  density  of  the  low-velocity  layer  on  the  flexural  wave 
group  velocity.  Three  near-overlapping  curves.  Top;  p„  =  0.6  x  IQ’  kg/m’  , 
middle:  =  0.8  x  10'  kg/m’  ,  and  bottom;  1.0  x  10’  kg/m’  .  The  model 

parameters  are:  a  =  3500  m/s,  j3  =  1856  m/s,  H  =  2  m,  =  0.89  x  10’ 
kg/m’,  V  =  1005  m/s,  d  =  0.2  m,  c  =  1440  m/s,  p,„  =  1  x  10’  kg/m’. 
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Fig.  3.8  Effect  of  ice  shear  wave  speed  on  flexural  wave  group  velocity.  Top:  /3  = 
1856  m/s,  middle:  /3  =  1756  m/s,  and  bottom;  jS  =  1656  m/s.  The  model 
parameters  are:  a  =  3500  m/s,  H  =  2  m,  p,  =  0.89  x  10’  kg/m’,  v  =  1005 
m/s,  d  =  0.2  m,  p„  =  0.8  x  10’  kg/m’,  c  =  1440  m/s,  =  lx  10’  kg/m’. 
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SECTION  4 


4  EFFECT  OF  VISCOUS  WAVES  FROM  SUSPENDED  PARTICLES  UNDER  SEA 

ICE  ON  LOW-FREQUENCY  ACOUSTIC  WAVE  REFLECTIVITY 

The  Arctic  ice  cover  affects  the  propagation  of  underwater  acoustic  waves,  and 
contributes  to  the  characteristics  of  underwater  ambient  noise.  It  is  essential  to  accurately 
model  sea  ice  to  be  able  to  predict  long  range  underwater  acoustic  wave  propagation. 

Hunkins  [4.1]  revealed  discrepancies  between  Arctic  field  data  and  theoretical  flexural 
wave  velocity  dispersion  curves  from  sea  ice.  The  calculated  ice  thickness  from  flexural  wave 
dispersion  was  about  20  %  smaller  than  the  actual  thickness.  Also,  the  measured  sea  ice 
flexural  wave  group  velocity  for  frequencies  above  about  100  Hz  was  about  up  to  40  % 
smaller  than  theoretically  predicted.  In  the  previous  section  of  this  report  (Section  3),  the 
effect  of  a  low-velocity  layer  underneath  floating  sea  ice  on  flexural  wave  dispersion  was 
discussed.  Theoretical  and  experimental  results  were  presented  giving  an  explanation  for 
observed  discrepancies  between  Hunkins’  flexural  wave  field  data  [4.1]  and  theoretical 
predictions.  By  introducing  a  low-velocity  fluid  layer  underneath  floating  sea  ice,  theoretical 
results  were  obtained  matching  the  flexural  wave  dispersion  field  results  [4.2].  Regarding 
under  ice  acoustic  wave  reflectivity,  Yang  and  Votaw  [4.3]  obtained  Arctic  field  results 
revealing  frequency  dependent  under  ice  reflectivity  data  that  could  not  be  matched  with  any 
known  theoretical  modeling  results. 

Fig.  4.1  compares  Yang  and  Votaw’s  field  measurements  with  theoretical  reflectivity 
calculations  obtained  by  the  author  using  a  thin  low-velocity  fluid  layer  between  the  ice  and 
the  water  half-space  as  shown.  The  theoretical  results  were  obtained  using  the  model 
parameters  given  in  Fig.  4.1.  One  notices  that  the  theoretical  and  field  reflectivity  results  do 
not  coincide.  The  introduction  of  a  low-velocity  fluid  layer  under  the  ice  does  not  fully  satisfy 
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Fig.  4.1  Under  ice  reflectivity  field  data  compared  with  theoretical  calculations  from 

ice  model  with  low  velocity  fluid  layer  under  floating  plate.  0  =  22.5°. 


Sonoquest  JRC'36-94 


J.  R.  Chamuel 


21 


the  reflectivity  data.  The  equations  used  in  this  section  to  calculate  the  reflection  coefficient 
were  derived  by  Folds  and  Loggins  (4.4).  The  reflectivity  results  presented  in  Fig.  4.1  were 
determined  using  a  shear  wave  speed  /3,  =  1  m/s  instead  of  zero  in  the  low-velocity  layer  to 
be  able  to  compare  these  results  with  the  ones  described  below  (identical  reflectivity  plots  were 
obtained  with  /3,  =  0  or  1  m/s). 

McCammon  and  McDaniel  [4.5]  computed  frequency  dependent  under  ice  theoretical 
reflectivity  results  from  a  layered  ice  model.  Yew  and  Weng  [4.6]  modeled  sea  ice  as  a 
transversely  isotropic  porous  medium  and  included  a  thin  porous  skeletal  layer  with  large  pore 
size  and  higher  porosity.  Their  theoretical  results  are  compared  with  Yang  and  Votaw’s  field 
data  (Fig.  4.2).  Qualitative  agreement  is  observed  between  the  field  and  calculated  results. 
Stanton  et  al.  |4.7j  used  reflections  of  100-800  kHz  normally  incident  underwater  acoustic 
waves  to  characterize  the  permeable  dendritic  structure  near  the  ice/water  interface.  Williams 
et  al.  [4.8]  conducted  in-situ  experiments  to  determine  the  acoustic  properties  of  the  skeletal 
layer  at  the  water/ice  interface. 

Chin-Bing  [4.9],  Posey  et  al.  [4.10],  and  Lawrence  et  al.  [4.11]  calculated  high 
frequency  under  ice  reflectivity  from  a  layered  ice  model  with  a  water-to-ice  transition  layer. 
The  ice  model  parameters  used  by  Chin-Bing  are  given  in  Table  4.1.  Notice  the  low  shear 
wave  velocity  in  this  transition  layer  model  has  a  corresponding  negligible  shear  wave 
attenuation  coefficient  (Table  4.1). 


HO.  3.  Idealized  ic»'bc  slmcture  (Schwarz  and  Weeks,  1977). 

C.  H.  Y«w  and  X.  Wang:  Raflacdon  and  rafracdon  from  sea  foa  343 


HG.  17.  Comparison  with  experimenuJ  resuhs.  a  »  66.6  de^  (grazing  an¬ 
gle  «  23.4  deg),  elastic  ice  layer  {N  ^  \  m,  0).  porous  ice  layer 
(//-  I  0.1).  skeleton  layer  (A  »  »  0.5). 


Fig.  4.2  Underice  reflectivity  field  data  compared  with  theoretical  calculations  by  Yew 
and  Weng’s  [4.6]  based  on  sea  ice  modeled  as  a  transversely  isotropic  porous 
medium.  0  =  23.4°. 
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Table  4.1  Layered  ice  and  transition  layer  parameters  used  by  Chin-Bing 
(4.9).  Notice  negligible  shear  wave  attenuation  in  transition 
layer.  (Table  copied  from  Reference  4.9) 
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Lindrose  (4.12]  used  shear  ultrasonic  waves  to  monitor  the  curing  process  of  epoxy 
resin.  His  results  showed  that  before  curing,  the  epoxy  shear  wave  speed  was  very  low, 
however,  the  corresponding  shear  wave  attenuation  was  very  high  (Fig.  4.3).  As  the  epoxy 
cured,  the  shear  wave  speed  increased  and  the  shear  wave  attenuation  decreased  rapidly.  The 
author  believes  that  there  is  a  similarity  between  curing  epoxy  and  solidifying  ice.  The  author 
was  able  to  match  Yang  and  Votaw’s  reflectivity  results  (4.3]  by  incorporating  the  effect  of 
shear  viscous  waves  in  the  low-velocity  layer  under  the  ice  with  an  extremely  high  shear  wave 
attenuation  coefficient.  Figure  4.4  compares  Yang  and  Votaw’s  results  with  Chamuel’s 
numerical  results  (4.2]  calculated  with  grazing  angle  ©  =  22.5°.  Very  good  agreement 
between  the  predicted  and  measured  results  was  obtained  when  the  shear  wave  speed  of  1  m/s 
in  the  low-velocity  layer  was  replaced  with  an  assumed  complex  shear  wave  velocity  equal  to 
1  -  650i  m/s.  No  attempt  was  made  here  to  model  the  viscoelastic  properties  of  ice  particles 
suspensions.  The  rheology  of  suspensions  is  beyond  the  scope  of  the  present  effort. 


Fig.  5— Shear-wave  speed  and  relative  attenuation 
vs.  epoxy  cure  time ;  •  present  data;  +  preliminary 
data 


Fig.  4.3 


Shear  wave  speed  and  relative  attenuation  plotted  vs  cure  time 
(copied  from  Lindrose,  Experimental  Mechanics,  1978]. 
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THEORETICAL  - 

(ChamueU 

EXPERIMENTAL  ARCTIC  FIELD  DATA  — » 
(YanQ  &  Votaw,  JASA  1981) 


0  200  400  600  SOO  1000 


FREQUENCY  (Hz) 


Fig.  4.4  Under  ice  reflectivity  field  data  compared  with  theoretical  calculations  based 
on  ice  model  with  viscous  waves  in  layer  under  floating  ice.  6  =  22.5°. 


Yang  and  Votaw’s  reflectivity  plots  depend  on  frequency  and  angle  of  incidence  of  the 
acoustic  waves.  If  the  model  proposed  by  Chamuel  is  correct,  it  should  satisfy  the  reflectivity 
dependency  on  both  frequency  and  grazing  angle.  Three  reflectivity  curves  plotted  as  function 
of  frequency  with  three  grazing  angles  0  =  18.3°,  22.5°,  and  23.4°  from  Yang  and  Votaw’s 
results  are  compared  in  Fig.  4.5  with  theoretical  curves  obtained  with  the  assumed  parameters 
shown  in  Fig.  4.6.  The  thickness  of  the  low-velocity  layer  d  was  9.61  cm  for  both  the  23.4° 
and  22.5°  grazing  angle  cases,  d  was  increased  to  16.2  cm  to  generate  the  reflectivity  curve 
at  0  =  18.3°.  The  geometrical  configuration  of  the  field  experiment  conducted  by  Yang  and 
Votaw  caused  the  acoustic  waves  to  reflect  from  different  horizontal  locations  along  the 
ice/water  interface  as  the  grazing  angle  was  changed.  For  example,  as  0  was  changed  from 
18.6°  to  23.4°,  the  reflection  point  shifted  by  52  m  horizontally  along  the  water/ice  interface. 
In  the  Arctic,  the  thickness  of  the  low-velocity  layer  is  range  dependent.  It  is  reasonable  to 
assume  that  the  thickness  of  the  low-velocity  layer  with  viscous  waves  changed  in  thickness 
from  9.6  cm  to  16.2  cm  between  the  two  locations.  Considering  that  there  are  several 
variables  and  unknowns  in  this  problem,  it  is  remarkable  to  have  been  able  to  match  the 
calculated  results  to  Yang  and  Votaw’s  field  results  for  the  first  time  both  in  frequency  and 
grazing  angle  (Fig.  4.5). 

The  presence  of  a  low-velocity  layer  supporting  shear  viscous  waves  under  sea-ice  may 
explain  the  observed  excessive  low-frequency  long-range  transmission  loss.  Future  Arctic 
field  experiments  need  to  include  in-situ  measurements  of  viscous  waves  in  a  thin  layer  near 
the  water/ice  interface. 
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Fig.  4.5  Underice  reflectivity  field  data  [4.3]  compared  with  theoretical  results  based  on 
ice  model  with  assumed  viscous  layer  under  floating  ice. 


Fig.  4.6  Ice  model  parameters. 
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SECTION  5 


5  DISPERSION  OF  ANTISYMMETRIC  EDGE  WAVES  ALONG  APEX  OF 
TRUNCATED  ELASTIC  WEDGE 


5. 1  Introduction 

Edge  waves  propagating  along  the  apex  of  an  elastic  wedge  were  discovered  in  1972 
independently  by  Lagasse  [5.1j  and  by  Maradudin  et  al.  15.2)  .  These  edge  waves  will  be 
referred  to  as  wedge  waves  in  this  report  in  order  not  to  confuse  them  with  flexural  edge 
waves  propagating  along  the  edge  of  a  plate  |5.3-5.5].  The  results  presented  in  this  section 
are  limited  to  the  study  of  antisymmetric  wedge  waves.  Fundamental  understanding  of  wedge 
waves  is  essential  for  nondestructive  testing,  underwater  acoustics,  acoustoelectronics, 
seismology,  acoustic  imaging,  and  noise  control.  Wedge  waves  can  propagate  for  example 
along  edges  of  propeller  blades,  finite  crystals,  cutting  tools,  mountains,  Arctic  ice,  cracks, 
airfoils,  underwater  fins,  and  ridges. 

Lagasse  predicted  the  existence  of  wedge  waves  using  the  finite  element  method  [5.6). 
Fig.  5.1  shows  Lagasse’s  results  revealing  the  computer-generated  wave  motion  along  the 
wedge  apex  and  the  velocity  of  antisymmetric  wedge  modes  plotted  as  function  of  wedge  apex 
angle  0.  Lagasse  obtained  an  empirical  relationship  for  the  antisymmetric  wedge  mode 
velocity 


V  =  Vk  sin  (n0),  n  =  1,  2,  3,  ...,  n0  <  90° 

where  is  the  Rayleigh  wave  velocity  on  a  half-space,  and  n  is  an  integer  representing  the 
mode  order.  The  antisymmetric  wedge  modes  are  nondispersive  and  the  number  of  existing 


the  three  lowest  order  nondispersive  wedge  waves.  Curves  were  com¬ 
puted  using  fourth-order  finite  elements  for  Duralumin  175.  Circles 
represent  experimental  results. 


Fig.  5.1  Wave  motion  and  velocity  of  antisymmetric  flexural  wedge  waves  plotted  as 
function  of  apex  angle  0.  [Copied  from  P.E.  Lagasse,  Reference  5.6). 
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modes  depends  on  0.  As  the  wedge  apex  angle  decreases,  the  number  of  antisymmetric  wedge 
modes  increases.  Experimental  results  using  broadband  pulses  were  presented  in  References 
|5.5  and  5.13].  Quantitative  studies  |5.4,  5.7,  5.8|  revealed  that  a  more  accurate  wedge  mode 
velocity  was  obtained  when  in  the  above  equation  was  replaced  with  the  edge  Rayleigh 
wave  velocity  V^r-;  existing  along  the  edge  of  a  quarter-space  (0  =  90°). 

Other  studies  on  wedge  waves  have  been  reported  in  several  papers.  Krylov  and 
Raguzina  15.9]  investigated  the  scattering  of  wedge  waves  by  notches.  Krylov  and  Shanin 
15.10]  determined  the  effect  of  anisotropy  on  wedge  waves.  Krylov  (5.11)  investigated  curved 
wedges  and  wedges  with  arbitrary  cross-section.  Krylov  and  Parker  [5.12]  developed  a  theory 
for  nonlinear  interactions  of  wedge  waves.  Chamuel  [5.13]  investigated  the  effect  of  water 
loading  on  wedge  waves,  range-dependent  apex  angle  [5. 7-5. 8]  and  immersed  90°  wedge 
[5.5],  Jia  et  al.  [5.14  -  5.15]  observed  unexplained  dispersive  behavior  of  flexural  wedge 
waves  using  laser  generation  and  detection. 

McKenna  et  al.  [5.16]  and  Lagasse  [5.17]  obtained  theoretical  and  numerical  results 
demonstrating  the  influence  of  apex  truncation  on  dispersion  of  antisymmetric  wedge  modes. 
Fig.  5.2  shows  a  normalized  phase  velocity  curve,  copied  from  Reference  [5.16],  for  the 
lowest  order  mode  on  a  truncated  wedge  plotted  as  function  of  X„/J,  where  is  the  distance 
from  the  "ideal"  apex  to  the  truncated  apex  surface,  and  /3  is  the  wavenumber.  Quantitative 


Fig.  6.  Normalized  phshe  velocity  for  the  lowest  order  mode  on  a 
truncated  wedge  as  functiona  of  XtS  in  the  interval  (0,  1 1  for  ►  -  0  0 
0.1,  0.17,  0.25,  0.333,  0.4,  0.6. 


Fig.  5.2  Phase  velocity  for  the  lowest  order  anitymmetric  flexural  wedge  mode 
propagating  along  apex  of  truncated  wedge.  [Copied  from  McKenna  et  al.. 
Reference  5.16]. 
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experimental  results  obtained  using  laser  ultrasonics  [5. 14-5. 15]  produced  dispersion  behavior 
of  antisymmetric  wedge  modes  not  predicted  by  theoretical  models.  The  purpose  of  this 
section  is  to  present  quantitative  experimental  results  characterizing  the  dispersion  of 
antisymmetric  wedge  waves  propagating  along  wedges  with  truncated  apex  using  contactless 
electromagnetic  acoustic  detection.  Contactless  electromagnetic  acoustic  detection  of 
antisymmetric  wedge  waves  provided  dispersion  results  in  agreement  with  the  theoretical 
predictions  by  McKenna  et  al.(5.16]. 

5.2  Experimental  setup 

Aluminum  6061  was  used  to  fabricate  a  wedge  with  an  apex  angle  0  =  26.7°.  The 
dimensions  of  the  wedge  were  chosen  to  separate  in  time  any  wave  reflections  from  the 
surfaces  and  corners  of  the  wedge  base  and  sides.  Ultrasonic  wedge  waves  were  generated 
with  a  1.5  mm  diameter  compressional  piezoelectric  transducer  positioned  to  excite  flexural 
antisymmetric  wedge  modes.  A  contactless  broadband  electromagnetic  acoustic  transducer  was 
used  to  detect  the  wedge  waves  propagating  along  the  apex.  The  principle  of  electromagnetic 
acoustic  transducers  was  described  in  a  paper  by  Thompson  [5.19].  The  geometrical 
configuration  of  the  detector  was  a  U  shaped  loop  similar  to  the  one  described  in  Reference 
[5.20]  for  magnetostrictive  ultrasonic  models.  The  received  signal  was  displayed  on  a  digital 
oscilloscope.  Waveform  and  dispersion  data  were  obtained  from  the  initially  sharp  wedge 
(3  urn  apex)  representing  the  "ideal"  wedge.  The  apex  was  then  truncated  successively  in 
steps,  and  corresponding  dispersion  results  were  recorded. 

5.3  Results 

Waveforms  of  the  antisymmetric  flexural  wedge  modes  were  obtained  with  the  source 
excited  with  a  single  broadband  pulse.  Phase  velocity  dispersion  plots  were  obtained  from 
gated  sinusoidal  signals  (30-500  KHz)  and  wavelength  measurements  were  averaged  over  10-20 


Fig.  5.3  Experimental  setup  configuration  and  waveform  of  nondispersive  antisymmetric 

flexural  wedge  waves  from  sharp  aluminum  wedge  with  h,,  =  3  urn. 
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wavelengths.  The  received  waveform  from  the  sharp  wedge  (h^  =  3  um)  is  displayed  in  Fig. 
5.3.  Three  practically  nondispersed  pulses  were  received  representing  the  three  allowed  modes 
for  the  26.7“  wedge.  The  distance  X  between  source  and  receiver  was  10  cm.  Fig.  5.4  shows 
dispersed  waveforms  obtained  from  the  truncated  wedge  with  h„  =  88  um  at  X  =  14.5  and 
X  =  20.5  cm.  A  highly  dispersive  wavetrain  was  produced  by  the  truncated  wedge  with  h„ 
=  200  um  (Fig.  5.5).  The  single  broadband  pulse  was  transformed  into  this  highly  dispersive 
wavetrain. 


Fig.  5.4  Waveform  of  dispersed  slowest  antisymmetric  flexural  wedge  wave  from 
truncated  apex  (h^  =  88  um). 


Fig.  5.5  Waveform  of  highly  dispersed  antisymmetric  flexural  wedge  wave  from 
truncated  apex  (h^  =  200  um). 


Sonoquest  JRC’36-94 


31 


J.  R.  Chamuel 


Fig.  5. 6-5. 9  show  the  measured  dispersion  of  antisymmetric  wedge  waves  obtained 
from  the  26.7°  aluminum  wedge  truncated  with  h„  =  30,  88,  140,  and  200  um.  The  ratio  of 
truncated  wedge  phase  velocity  V  to  the  "ideal"  wedge  was  plotted  as  function  of 
frequency.  The  velocity  of  the  high  frequency  components  increased  substantially  as  h^  was 
Increased.  The  repeatability  of  the  velocity  data  points  was  better  than  10'4  up  to  about  500 
kHz.  Above  500  kHz,  the  signal-to-noise  ratio  dropped  drastically  and  the  accuracy  of  the 
data  points  obtained  from  the  truncated  wedge  with  h,,  =  30  um  was  questionable  (Fig.  5.6). 
Notice  the  vertical  scale  is  different  in  Fig.  5. 6-5. 9. 

In  order  to  compare  the  experimental  results  on  the  effect  of  apex  truncation  on  the 
dispersion  of  the  slowest  antisymmetric  flexural  wedge  mode  with  the  theoretical  predictions 
by  McKenna  et  al.  (5.16],  Fig.  5.10  was  generated.  The  ratio  V/V^  was  plotted  vs  Xj5  as 
presented  in  Reference  [5.16].  The  experimental  data  gathered  with  the  contactless 
electromagnetic  acoustic  transducer  compared  somewhat  favorably  with  the  theoretical  results. 
The  anomalous  laser  ultrasonics  experimental  data  described  in  Reference  [5.15]  with  0  =  30° 
are  replotted  in  Fig.  5.10. 


5.3  Conclusions 

Accurate  quantitative  results  were  presented  characterizing  the  dispersive  behavior  of 
antisymmetric  flexural  wedge  waves  propagating  along  the  apex  of  a  truncated  elastic  wedge. 
Contactless  electromagnetic  acoustic  detection  of  antisymmetric  wedge  waves  produced 
dispersion  results  generally  in  agreement  with  theoretical  predictions  by  McKenna  et  al.  [5.16]. 


h  =  30  um 

O 


Fig.  5.6  Measured  dispersion  of  antisymmetric  wedge  waves  obtained  from  truncated 
aluminum  wedge  (h,,  =  30  um). 
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h  =  88  um 
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Fig.  5.7  Measured  dispersion  of  antisymmetric  wedge  waves  obtained  from  truncated 
aluminum  wedge  (h„  =  88  um). 


140  um 


Fig.  5.8  Measured  dispersion  of  antisymmetric  wedge  waves  obtained  from  truncated 
aluminum  wedge  (h^  =  140  um). 
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Fig.  5.9  Measured  dispersion  of  antisymmetric  wedge  waves  obtained  from  truncated 
aluminum  wedge  (h„  =  200  um). 
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Fig.  5.10  Experimental  results  on  effect  of  apex  truncation  on  dispersion  of  the  slowest 
antisymmetric  flexural  wedge  mode  compared  with  theoretical  results  by 
McKenna  et  al.  [5.16].  The  anomalous  laser  ultrasonics  experimental  data 
with  6  =  30°  from  Reference  [5.15]  are  also  plotted. 
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SECTION  6 


6  ELASTIC  WEDGE  WITH  RANGE-DEPENDENT  APEX  ANGLE  AND  CROSS 
SECTION 

The  following  two  papers  summarize  new  findings  obtained  on  edge  wave  propagation 
along  wedges  with  range-dependent  apex  angle  and  cross  section. 

RANGE-DEPENDENT  ANTISYMMETRIC  EDGE  WAVES  ALONG  CURVED 
ELASTIC  WEDGE  WITH  VARIABLE  CROSS  SECTION 

Jacques  R.  Chamuel 

Sonoquesl  Advanced  Ultrasonics  Research 

P.O.  Box  153 

Wellesley  Hills,  MA  02181-5339 
INTRODUCTION 

Flexural  waves  in  isotropic  plates  are  known  to  be  dispersive  [1].  Edge  flexural 
waves  propagating  along  the  edge  of  a  semi-infinite  plate  are  slower  than  conventional 
flexural  waves  and  are  also  dispersive  [2],  The  only  known  nondispersivc  flexural 
waves  are  antisymmetric  edge  waves  propagating  along  the  apex  of  a  sharp  elastic  wedge 
discovered  by  Lagasse  [3].  The  present  paper  focuses  on  antisymmetric  edge  waves 
propagating  along  the  apex  of  an  elastic  wedge  with  a  range-dependent  apex  angle. 

Edge  waves  exist  in  a  large  number  of  geometrical  configurations  and  may  propagate 
along  ridges,  faults,  propeller  blades,  sea-ice  features,  mountain  apex,  oblique  cracks  in 
cores,  boreholes  intersecting  oblique  boundaries,  ripple  marks  in  sediments,  and  helical 
drills.  Edge  waves  can  be  generated  by  conventional  acoustic  sources  or  by  acoustic 
emission  from  fractures.  This  research  is  aimed  at  developing  fundamental  physical 
understanding  of  coupling  between  underwater  acoustic  waves  and  elastic  waves  in  finite 
solids  with  oblique  boundaries  and  range-dependent  effective  elastic  properties.  Results 
on  edge  waves  along  immersed  wedges  were  presented  elsewhere  [4-5]. 

WEDGE  EDGE  WAVES 

Along  the  apex  of  an  clastic  wedge,  nondispersivc  localized  antisymmetric  edge 
waves  (also  known  as  line  acoustical  waves  [6])  exist.  Their  particle  displacement 
decays  exponentially  away  from  the  boundary.  Of  all  the  different  types  of  elastic  waves 
known,  nondispersivc  wedge  edge  flexural  waves  have  the  highest  energy  density  and 
largest  strain.  Lagasse  [7]  introduced  a  simple  empirical  equation  for  the  phase  velocity 
V  of  each  antisymmetric  edge  mode 

V  *  Vr  sin  (m0),  m  =  1,2,3,...,  m0  <  90®  (1) 

where  is  the  Rayleigh  wave  velocity,  6  is  the  wedge  apex  angle,  and  m  is  the  mode 
order.  The  author  found  that  a  more  accurate  expression  is  obtained  when  in 
equation  (1)  is  replaced  with  the  slower  edge  Rayleigh  wave  velocity  existing  along 
the  edge  of  a  solid  quarter  space  [8-9],  V  =  sin  (mS). 
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Fig.  1.  (a)  Sketch  of  cylindrical  model,  (b) 
Measured  edge  angle  along  elliptical  edge. 


EXPERIMENTAL  SETUP 

The  experiments  were  performed  on  two  basic  geometrical  configurations:  a) 
truncated  obliquely  cut  cylinder  and  (b)  obliquely  truncated  borehole.  The  models  were 
machined  from  6061  aluminum.  The  diameter  D  of  the  models  was  10. 1 19  cm  and  the 
cut  was  at  an  angle  of  (i>  =  22.5°.  Broadband  pulsed  ultrasonic  edge  waves  were 
generated  and  detected  using  1.5  mm  piezoelectric  compressional  wave  transducers. 

RESULTS  AND  DISCUSSION 

Fig.l  shows' a  sketch  of  the  cylindrical  model  and  the  measured  edge  angle  6  as 
function  of  arc  length  S.  The  reference  point  S  =  0  was  located  at  $  =  90°  with  the 
positive  direction  towards  the  acute  apex  as  shown.  0  continuously  varied  between  22.5° 
and  157.5°.  The  velocity  of  antisymmetric  edge  modes  were  calculated  as  function  of 
apex  angle  for  the  range-independent  flat  wedge  case  (Fig.  2(a)).  The  ratio  of  V/Vre 
was  plotted  as  function  of  0  up  to  0  =  90°.  Note  that  the  number  of  wedge  modes 
increased  as  6  decreased.  The  waveform  of  edge  waves  propagating  along  the  direction 
of  increasing  6  were  compared  with  waves  along  decreasing  6  (Fig.  2(b)).  The  source 
was  located  at  S  =  0  (0  =  90°).  At  S  =  -15  cm,  6  was  157.5°  and  only  one  edge 
wave  mode  propagating  near  the  Rayleigh  wave  velocity  was  observed  (Fig.  2(b)  bottom 
trace).  At  S  =  +  15  cm,  6  was  22.5°  and  three  modes  were  detected.  During  the 
mode  conversion  process,  most  of  the  energy  was  concentrated  in  the  slowest  mode. 


Fig.  2.  (a)  Velocity  of  antisymmetric  wedge  modes  as  function  of  flat  wedge  apex  angle, 
(b)  Comparison  of  edge  waves  at  S  =  -H5  cm  (top  trace)  and  -  15  cm  (bottom  trace). 
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Fig.  3.  Measured  edge  angle  vs  arc  length  compared  with  approximated  linear  segments. 


A  plot  of  the  measured  edge  angle  as  function  of  arc  length  is  shown  in  Fig.  3. 
Two  linear  segments  were  used  to  approximate  5(S).  For  the  first  quarter  of  the  ellipse^ 

e  =  -7.81  IS  +  t/2  ,  (2) 


Taking  x  as  the  principal  axis  of  the  ellipse  and  x  =  0  at  its  center,  the  x-iniercepts  are 
±  D  /  2  sin  .  The  edge  wave  velocity  was  calculated  using  the  approximate  linear 
segments  and  the  modified  Lagasse’s  equation  with  Vj<,:=  2886  m/s  (Fig.  4)  instead  of 
Vr=2951  m/s.  The  arc  length  was  determined  from  the  expression 


sin^(t)(- 


4x  ^sin^(J) 


(3) 


and  the  average  edge  wave  velocity  (Fig,  4)  was  computed  from 


(4) 


Fig.  4.  Calculated  edge  wave  velocity  and  average  velocity  for  the  slowest  mode. 

1425 


Sonoquest  JRC-36-9A 


39 


J.  R.  Chamuel 


EXPERIMENTAL  DATA  o  o  o  o 
APPROXIMATE  LINEAR  MODEL  - 


0 


0 


0.05  0.1  0.15  0.2 

ARC  LENGTH  S  (m) 


Fig.  5.  Comparison  of  measured  and  calculated  time  of  arrival  of  slowest  edge  wave 
component. 


Fig.  5  reveals  the  nonlinear  characteristics  of  edge  wave  propagation  time  and  compares 
the  measured  and  calculated  time  t  =  S  /  v  for  the  slowest  edge  wave  pulse.  Edge 
waves  are  elastic  waves  localized  near  the  edge  of  the  solid  and  decay  exponentially 
inside  the  solid.  In  these  range-dependent  models,  the  degree  of  localization  of  wedge 
waves  is  also  variable  since  it  depends  on  the  edge  wave  velocity.  Near  6  =  22.5®,  the 
edge  wave  is  more  localized  than  at  6  =  90®,  Fig.  6  demonstrates  the  decay  of  the  edge 
wave  with  distance  from  the  edge  by  displaying  the  spectrum  of  the  slowest  edge  wave 
component  as  function  of  the  distance  h  from  the  edge.  The  high  frequency  components 
were  attenuated  more  rapidly  than  the  low  frequency  components  as  h  was  increased. 


Fig.  6.  Decay  of  edge  wave  with  distance  from  edge  demonstrating  localization.  Top:  at 
6  =  90®,  Bottom:  at  0  =  22.5®.  Note  greater  localization  near  acute  edge  angle. 
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Fig,  7.  Waveforms  recorded  at  different  ranges  along  the  elliptical  path. 


A  collection  of  ultrasonic  waveforms  recorded  at  different  ranges  S  =  5,  10,  15, 
20,  25,  and  30  cm  is  presented  in  Fig.  7.  In  addition  to  the  edge  waves,  a  direct 
Rayleigh  wave  arriving  before  the  edge  w'aves  and  propagating  across  the  Hat  elliptical 
surface  was  detected.  The  waveforms  demonstrate  the  continuous  mode  conversion  as 
the  edge  waves  propagated  from  a  90°  wedge  angle  to  22.5°  and  again  to  90°.  The 
relative  amplitude  of  the  edge  wave  particle  velocity  in  the  direction  normal  to  the  flat 
plane  was  measured  and  plotted  as  function  of  S  (Fig.  8).  The  amplitude  first  increased 
and  then  dipped  near  6  =  12°-18°.  The  waveshape  of  the  slowest  edge  mode  near  the 
acute  apex  is  closely  examined  in  Fig.  9. 


1  ' 


0  to  2u  JO  40 


^  S  (cm) 

Fig.  8.  Measured  amplitude  of  the  slowest  wave  component  along  the  elliptical  edge. 
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Experimental  results  from  the  elliptical  edge  of  a  truncated  borehole  gave 
practically  the  same  type  of  edge  wave  characteristics  in  the  used  frequency  range. 

Time  of  arrival  data  from  the  borehole  overlapped  the  results  of  Fig.  3(c)  from  the 
truncated  cylinder.  The  direct  Rayleigh  wave  component  crossing  the  minor  axis  of  the 
ellipse  was  obviously  absent  in  the  borehole  case.  Fig.  10  compares  the  waveforms 
obtained  from  the  cut  solid  cylinder  and  oblique  borehole.  A  high  frequency  component 
was  observed  at  the  onset  of  the  waveform  from  the  cylinder.  Further  studies  on  edge 
waves  are  in  progress. 


Fig,  10.  Comparison  of  slowest  antisymmetric  edge  wave  from  oblique  borehole  and 
solid  cylinder  cut  at  22,5°.  Source  located  at  S  =?  0  and  receiver  at  at  S  =  30.16  cm. 
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Fig.  1 1.  Range-dependent  antisymmetric  edge  waves  along  propeller  blade  edge. 


A  preliminary  experiment  was  conducted  on  the  edge  of  a  model  airplane  plastic 
propeller  (Windsor  Propeller  No.  1080)  to  demonstrate  the  existence  of  edge  waves  (Fig. 

1 1(a).  Source  and  receiver  transducers  were  placed  along  an  edge  of  the  propeller  blade 
as  shown  separated  by  5  cm.  The  detected  antisymmetric  edge  wave  was  large  and 
dispersive.  Edge  wave  dispersion  was  due  to  truncation  of  the  propeller  edge. 

Dispersion  ettects  of  truncating  the  apex  of  a  wedge  were  previously  investigated  [10- 
1 1).  Localization  of  the  antisymmetric  edge  wave  was  demonstrated  by  the  large 
attenuation  of  the  wave  as  the  receiver  R2  was  placed  7  rnm  away  from  the  edge  (Fig. 
11(b)). 

SUMMARY  AND  CONCLUSIONS 

This  paper  presented  the  first  results  on  range-dependent  broadband  pulsed 
antisymmetric  edge  waves  propagating  along  the  sharp  edge  of  a  curved  elastic  wedge 
with  variable  cross  section.  Experimental  results  from  truncated  oblique  cylinders  and 
truncated  oblique  boreholes  revealed  that  the  velocity,  amplitude,  and  number  of  existing 
antisymmetric  edge  modes  continuously  changed  as  the  waves  propagated  along  the 
range  dependent  edge.  The  antisymmetric  edge  waves  remained  nearly  nondispersive 
over  a  wide  frequency  range.  For  a  10  cm  diameter  aluminum  cylinder  obliquely  cut  at 
22.3°,  the  edge  angle  around  the  elliptical  path  varied  almost  linearly  as  function  of  arc 
length  between  22.5°  and  90°.  As  the  apex  angle  increased,  the  number  of  edge  modes 
decreased.  An  example  was  given  on  edge  waves  propagating  along  the  edge  of  a 
propeller.  In  view  ot  space  limitations,  several  references  on  edge  w'aves  along  fixed 
angle  wedges  were  not  mentioned.  The  reader  may  refer  to  publications  by  Krylov  [12], 
The  new  findings  contribute  to  the  physical  understanding  of  3-D  range  dependent 
liquid/solid  seismoacousiic  wave  phenomena  and  propeller  noise. 
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EDGE  WAVES  ALONG  IMMERSED  ELASTIC  ELLIPTICAL 
WEDGE  WITH  RANGE  DEPENDENT  APEX  ANGLE 

Jacques  R,  Chamuel 

Sonoquest  Advanced  Ultrasonics  Research 
P.O.  Box  153,  Wellesley  Hills.  MA  02181-5339 


ABSTRACT 

Experimental  results  are  presented  revealing  the 
characteristics  of  wedge  waves  propagating  along  the  edge 
of  free  and  immersed  curved  elastic  wedges  with  variable 
cross  section  and  range-dependent  apex  angle.  Broadband 
pulsed  wedge  waves  from  the  elliptical  edge  of  obliquely 
cut  solid  cylinders  and  boreholes  arc  compared.  The 
velocity  and  number  of  existing  edge  modes  continuously 
vary  along  the  wave  path.  Leaky  edge  waves  radiate  in 
the  water  with  continuously  varying  angles  depending  on 
the  local  edge  wave  velocity.  Antisymmetric  edge  waves 
remained  nearly  nondispersive  over  a  wide  frequency 
range.  The  new  findings  provide  physical  insight  in 
understanding  range-dependent  liquid/solid  wave 
phenomena  and  propeller  noise. 

INTRODUCTION 

Nondispersive  flexural  wedge  waves  propagating  along  the 
apex  of  an  ideal  sharp  clastic  wedge  were  discovered  in 
1972  by  Lagasse  |1)  and  have  been  studied  by  several 
researchers  since  [2].  The  terms  edge  or  wedge  waves  will 
be  used  interchangeably  in  the  text.  The  particle 
displacement  *  of  these  localized  nondispersive 
antisynunetric  edge  waves  decays  exponentially  away 
from  the  boundary.  Of  all  the  different  types  of  clastic 
waves  known,  nondispersive  wedge  flexural  waves  have 
the  highest  energy  density  and  largest  strain.  The  present 
paper  focuses  on  curved  elastic  wedges  with  variable  cross 
section  and  range-dependent  apex  angle.  Wedge  waves 
exist  in  a  large  number  of  geometrical  configurations  and 
may  propagate  along  ridges,  sea-ice  features,  mountains, 
oblique  cracks  in  cylindrical  cores,  oblique  boundaries 
intersecting  boreholes,  ripple  marks  in  sediments,  helical 
drills,  and  propeller  blades.  Edge  waves  may  be 
generated  by  conventional  acoustic  sources  or  by  acoustic 
emission  from  cracking.  This  research  is  aimed  at 
developing  fundamental  physical  understanding  of  the 
coupling  between  underwater  acoustic  waves  and  elastic 


waves  in  finite  solids  with  oblique  boundaries  and  range- 
dependent  effective  elastic  properties.  Other  results  on 
wedge  waves  were  presented  elsewhere  (4-5). 

Lagasse  (1)  introduced  a  simple  empirical  equation  for  the 
phase  velocity  V  of  each  antisymmetric  edge  mode 

V  =  V,^  sin  (m5),  m  *=  1,2,3,...,  m^  <  90®  (1) 

where  is  the  Rayleigh  wave  velocity,  6  is  the  wedge 
apex  angle,  and  ra  is  the  mode  order.  According  to  H). 
the  velocity  of  each  antisymmetric  mode  is  independent 
of  frequency,  and  it  decreases  as  the  wedge  angle  is 
decreased.  The  author  found  that  a  more  accurate 
expression  is  obtained  when  V,j  in  equation  (1)  is 
replaced  with  the  slower  edge  Rayleigh  wave  velocity 
Vn£  existing  along  the  edge  of  a  solid  quarter  space  [5], 

V  =  sin  (m^). 

In  addition  to  the  antisymmetric  wedge  modes, 
symmetric  edge  waves  with  a  velocity  close  to  the 
Rayleigh  wave  can  also  exist. 


EXPERIMENTAL  SETUP 

The  experiments  were  performed  on  three  basic 
geometrical  configurations:  a)  a  reference  mode! 
consisting  of  a  conventional  triangular  wedge  with  a 
fixed  apex  angle,  b)  a  truncated  obliquely  cut  cylinder 
and  (c)  an  obliquely  truncated  borehole.  The  models 
were  machined  from  6061  aluminum.  The  triangular 
wedge  had  an  apex  angle  ©  »  26.72^.  The  borehole 
and  cylinder  had  the  same  diameter  D  *  10. 1 19  cm  and 
the  cut  was  at  an  angle  of  ^  =  22.5®.  Broadband 
pulsed  ultrasonic  edge  waves  were  generated  and 
detected  using  1.5  mm  diameter  piezoelectric 
compressional  wave  transducers. 
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RESULTS  AND  DISCUSSION 
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I.  Fixed  Wedge  Angle 

Antisymmetric  edge  waves  were  generated  in  the 
triangular  wedge  with  the  transducers  oriented 
perpendicular  to  the  apex  to  transversally  excite  the 
flexural  wedge  waves  Fig.  1(a).  Three  antisymmetric 
modes  were  detected,  however  the  slowest  mode  (mode 
1)  was  very  dispersive  mainly  due  to  the  finite 
truncation  of  the  wedge  apex  (30  /im).  The  effect  of 
apex  truncation  on  the  dispersion  of  flexural  wedge 
waves  was  previously  investigated  [7*8].  By  sharpening 
the  wedge,  the  dispersion  was  significantly  reduced 
(Fig.  1(b)).  The  triangular  wedge  was  fabricated  by 
diagonally  cutting  and  machining  a  7.5  cm  thick 
aluminum  plate.  The  measured  Rayleigh  wave  velocity 
on  one  face  of  the  wedge  was  0.56%  smaller  than  the 
other  side.  Observed  dispersive  behavior  of  laser 
generated  antisymmetric  wedge  waves  were  reported 
(9].  Synunetric  edge  waves  were  generated  with  the 
transducers  onenied  parallel  to  the  apex  (Fig.  2).  The 
first  observed  pulse  was  a  longitudinal  wave  followed  by 
a  large  symmetric  edge  Rayleigh  wave.  The 
source/receiver  distance  x  was  25.5  cm. 


d3us 


II.  RangeOependent  Wedge  Angle 

Fig. 3  shows  a  sketch  of  the  obliquely  cut  cylindrical 
model.  The  measured  edge  angle  6  as  function  of  arc 
length  S  is  plotted  in  Fig.  4.  6  varied  continuously 
between  22.5®  and  157,5®.  The  reference  point  S  =  0 
was  located  at  ^  =  90®  with  the  positive  direction 
towards  the  acute  apex  as  shown.  The  source  was 
located  at  S  =  0  on  the  edge.  The  receiver  was  placed 
on  the  flat  elliptical  surface  away  from  the  edge  at  a 
distance  r  -  9  cm.  A  small  direct  Rayleigh  wave  was 
observed  as  expected  (Fig.  5(a)).  The  receiver  was  then 
moved  to  S  =  9  cm  on  the  edge  to  detect  the  edge  wave 
displayed  in  Fig.  5(b)),  This  experiment  demonstrates 
the  relative  amplitude  of  the  direct  Rayleigh  wave  on  the 
flat  surface  and  the  antisymmetric  edge  wave. 


Fig.  1.  Effect  of  wedge  apex  truncation  on 

antisymmetric  edge  waves.  Top:  30  ^m  wedge  tip  Fig.  3.  Sketch  of  cylindrical  model, 

thickness.  Bottom:  sharp  wedge  tip. 
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Fig.  5.  Comparison  of  direct  R  wave  and  edge  wave. 


A  collection  of  ultrasonic  waveforms  recorded  at 
different  ranges  S  «  5.  10,  15  ,  20.  25,  and  30  cm  is 
presented  in  Fig.  6.  In  addition  to  the  edge  waves,  a 
direct  Rayleigh  wave  arriving  before  the  edge  waves  and 
propagating  across  the  flat  elliptical  surface  was 
detected.  The  waveforms  demonstrate  the  continuous 
mode  conversion  as  the  edge  waves  propagated  from  a 
90®  wedge  angle  to  22.5®  and  again  to  90®. 


Two  linear  segments  were  used  to  approximate  $(S). 
For  the  first  quarter  of  the  ellipse, 

e  «  -7.81  IS  +  t/2  . 


Taking  x  as  the  principal  axis  of  the  ellipse  and  x  =  0 
at  its  center,  the  x-intercepts  are  ±  D  /  2  sin  o  .  Tr.e 
edge  wave  velocity  was  calculated  using  the  approxirnate 
linear  segments  and  the  modified  Lagaisse’s  equation 
with  V^£=  2820  m/s  CFig.  7)  instead  of  Vj  =  2932  m,  s 
The  arc  length  was  determined  from  the  expression 


4A:^sin^(}) 


•1)  \dx  . 


Fig.  6.  Waveforms  recorded  at  different  ranges. 
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and  the  average  edge  wave  velocity  (Fig.  7)  was 
computed  from 
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0.0$  0.1  0.1$  0.3  0.3$  0.3 


Fig.  7.  Calculated  edge  wave  velocity  and  average 
velocity  for  the  slowest  mode. 

Fig.  8  reveals  the  nonlinear  characteristics  of  edge  wave 
propagation  time  and  compares  the  measured  and 
calculated  time  I  =  S  /  v  for  the  slowest  edge  wave 
pulse.  Fig.  9  displays  the  calculated  local  velocity  of 
the  three  antisymmetric  modes.  Edge  waves  are  elastic 
waves  localized  near  the  edge  of  the  solid  and  decay 
exponentially  inside  the  solid.  In  these  range*dependent 
models,  the  degree  of  localization  of  wedge  waves  is 
also  variable  since  it  depends  on  the  edge  wave  velocity. 
Near  G  =  22.5®,  the  edge  wave  is  more  localized  than 


Fig.  9.  Calculated  local  velocity  of  three  antisymmetric 
edge  modes. 

Experimental  results  from  the  elliptical  edge  of  a 
truncated  borehole  gave  practically  the  same  type  of 
edge  wave  characteristics  in  the  used  frequency  range. 
Time  of  arrival  data  from  the  slowest  borehole  mode 
nearly  overlapped  the  results  of  Fig.  8  from  the 
truncated  cylinder.  Fig.  10  compares  the  antisymmetric 
edge  waves  from  the  cylindrical  model  (Fig.  10(a))  and 
borehole  (Fig.  10(b))  with  S  =  15.08  cm  (one  quarter 
the  perimeter).  The  experiment  was  repeated  with  the 
receiver  at  S  =  30.16  cm  (Fig.  ll).  The  direct 
Rayleigh  wave  component  crossing  the  minor  axis  of  the 
ellipse  was  obviously  absent  in  the  borehole  case.  Also 
the  first-arriving  antisymmetric  edge  waves  were 
severely  attenuated  in  the  borehole  model.  Thjs  is  due 
to  leakage  of  these  waves  into  body  waves  due  to  the 
curvature  of  the  edge. 


Fig.  8.  Comparison  of  measured  and  calculated  time  of 
arrival  of  slowest  antisymmetric  edge  wave  component. 


Fig.  10.  Comparison  of  borehole  and  solid  cylinder 
waveforms  of  antisymmetric  edge  waves  (S  =  15  cn 
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Fig.  11.  Comparison  of  borehole  and  solid  cylinder 
waveforms  of  antisymmetric  edge  waves  (S  =  30  cm). 


Experiments  were  conducted  on  the  solid  cylinder  and 
borehole  models  to  demonstrate  the  effect  of  water 
loading  on  the  slowest  antisymmetric  edge  wave  (Fig 
12-13).  Water  loading  decreased  the  velocity  of  the 
slowest  flexural  mode  and  also  it  introduced  dispersion 
and  attenuation  as  shown.  In  these  controlled 

expenments,  the  water  level  was  increased  lo  wet  a 


range  was  fixed  30.16  cm.  The  arc  length  S' in  the 
figures  corresponds  lo  the  distance  between  the  source 
and  the  water  level. 


The  velocity  of  the  antisymmetric  edge  waves  vary  from 
point  to  point  around  the  elliptical  path.  The  supersonic 
antisymmetric  wave  components  leak  their  energy  to  the 
water  with  different  angles  that  depend  on  the  local 
velocity.  Fig.  14  shows  the  calculated  radiation  angle 
for  the  fastest  mode  as  function  of  arc  length  assuming 
the  free  edge  wave  velocity. 


Fig.  12.  Effect  of  water  loading  on  the  slowest 
antisymmetric  edge  wave  component  in  solid  cylindneal 
model. 
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Fig.  13.  Effect  of  water  loading  on  the  slowest 
antisynunetric  edge  wave  component  along  borehole. 


so 


0  0.02  0.04  0.06  0.08  0.1 


s 

Fig.  14.  Calculated  radiation  angle  of  supersonic 
antisynunetric  edge  waves  as  function  of  position  along 
the  elliptical  wavepath.  Aluminum/water  model  based  on 
free  model  velocity. 

SUMMARY  AND  CONCLUSIONS 

New  findings  were  presented  on  range-dependent 

Sonoquest  JRC-36-94 


broadband  pulsed  antisymmetric  edge  waves  propagating 
along  the  sharp  edge  of  a  curved  elastic  wedge  with 
variable  cross  section.  Experimental  results  from 
truncated  oblique  solid  cylinders  and  truncated  oblique 
boreholes  revealed  that  the  velocity,  amplitude,  and 
number  of  existing  antisymmetric  edge  modes 
continuously  changed  as  the  waves  propagated  along  the 
range  dependent  edge.  The  antisymmetric  edge  waves 
remained  nearly  nondispersive  over  a  wide  frequency 
range.  For  a  10  cm  diameter  aluminum  cylinder 
obliquely  cut  at  22.5®,  the  edge  angle  around  the 
elliptical  path  varied  almost  linearly  as  function  of  arc 
length  between  22.5*  and  90®.  As  the  apex  angle 
increased,  the  number  of  edge  modes  decreased.  The 
effect  of  water  loading  on  antisymmetric  edge  waves 
was  demonstrated.  The  new  findings  contribute  to  the 
physical  understanding  of  3-D  range  dependent 
iiquid/solid  seismoacoustic  wave  phenomena  and 
propeller  noise. 
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SECTION  7 


7  SHEAR  WAVE  COUPLING  TO  ICE  CORES 

One  of  the  greatest  difficulties  encountered  in  experimental  studies  of  sea  ice  is  shear 

wave  coupling  between  transducer  and  ice  17. 1|.  Efficient  coupling  of  ultrasonic  waves  can 
be  achieved  by  bonding  the  transducer  to  the  surface  of  the  solid,  however,  bonding  may  not 
be  practical  in  many  situations.  Coupling  of  compressional  waves  can  be  achieved  by  using 
a  liquid  couplant  such  as  water  between  the  transducer  and  ice.  Shear  waves  on  the  other  hand 
are  not  transmitted  through  water  and  a  viscous  liquid  such  as  honey  for  example  may  be  used 
to  couple  shear  waves  to  ice.  Commercially  available  shear  wave  couplants  (e.g.  Panametrics 
SWC  couplant)  are  not  practical  for  field  testing  Arctic  ice  cores  at  very  low  temperatures 
because  of  their  high  viscosity  and  because  the  ice  samples  become  contaminated.  Dry 
coupling  between  the  ultrasonic  transducer  and  the  solid  is  desirable  in  several  applications 
such  as  testing  ice  cores,  however,  dry  coupling  between  conventional  ultrasonic  shear  wave 
transducers  and  ice  allows  little  energy  to  be  transmitted  to  the  ice. 

The  author  introduced  a  new  technique  for  dry  coupling  ultrasonic  shear  waves  between 
a  shear  transducer  and  ice  |7.21.  The  technique  generates  and  detects  shear  ultrasonic  waves 
using  an  abrasive  interface  layer  composed  of  fine  hard  particles  forming  a  sandpaper-like 
surface  that  causes  gripping  between  the  transducer  and  the  solid  as  a  small  force  is  applied 
to  the  transducer  normal  to  the  surface  of  the  solid. 

The  proposed  method  can  be  implemented  by  bonding  high  friction  hard  particles  (e.g. 
silicon  carbide,  carborundum,  diamond,  alumina  zirconia,  aluminum  oxide,  etc...)  directly  to 
the  smooth  surface  of  a  conventional  ultrasonic  shear  transducer,  or  to  a  wear  plate  on  the 
surface  of  the  transducer.  Another  approach  is  to  use  a  removable  double-sided  thin  sheet 
coated  on  both  sides  with  high  friction  hard  particles  to  be  inserted  between  a  conventional 
ultrasonic  shear  wave  transducer  and  the  solid. 

Experimental  results  verifying  the  proposed  technique  are  shown  in  Fig.  7.1.  Two 
identical  shear  wave  transducers  (Panametrics  V152,  1  MHz,  1"  diameter)  were  placed  on 
opposite  ends  of  a  10  cm  thick  ice  block  as  shown.  The  transducers  were  oriented  such  that 
their  polarization  was  in  the  vertical  direction  to  generate  and  detect  a  shear  vertical  wave  in 
the  ice  block.  The  source  was  excited  with  a  single  broadband  pulse  2  us  wide.  Fig.  7.1(a) 
shows  the  waveform  of  the  received  signal  (horizontal  time  scale  10  us/div.).  The  first- 
arriving  pulse  was  a  compressional  wave  and  no  significant  shear  wave  was  detected.  Silicon 
carbide  powder  grit  220  was  then  bonded  to  the  surface  of  two  transducers,  and  the  above 
experiment  was  repeated.  The  received  waveform  from  the  transducers  coated  with  silicon 
carbide  powder  is  displayed  in  Fig.  7.1(b).  The  first  arriving  pulse  is  a  compressional  wave 
and  the  second  large  pulse  is  the  desired  shear  wave. 
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compressional  shear  wave 

wave 


Fig.  7.1  Waveform  of  elastic  waves  generated  and  detected  in  ice.  (a)  using 
conventional  shear  transducer  with  smooth  surface,  (b)  using  shear  transducer 
coated  with  silicon  carbide  powder. 
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SECTION  8 


8  ULTRASONIC  CHARACTERIZATION  OF  ARCTIC  ICE  CORES 

In  order  to  develop  comprehensive  understanding  of  the  physical  properties  of  sea  ice, 
R.  Obrochta/ONR  proposed  the  use  of  a  portable  ultrasonic  system  for  field  measurements  of 
compressional  and  shear  wave  velocities  continuously  along  ice  cores  to  complement  other 
Arctic  ice  studies.  The  author  investigated  a  related  problem  on  variations  in  shear  wave 
velocity  in  cylindrical  ice  cores  across  two  orthogonal  diametrical  paths  at  a  given  ice  core 
depth.  The  purpose  of  this  section  is  to  present  laboratory  ultrasonic  results  revealing 
variations  in  shear  wave  velocity  across  two  orthogonal  diametrical  paths  in  ice  core  samples 
18.1-8. 2j. 


Sea-ice  core  samples  were  provided  by  W.  B.  Tucker/  U.  S.  Army  Cold  Regions 
Research  and  Engineering  Laboratory.  Shear  source  and  receiver  transducers  were  positioned 
opposite  each  other  across  the  core  diameter.  The  source  was  excited  with  a  2  us  broadband 
pulse  and  the  received  waveform  was  captured  by  a  digital  oscilloscope.  One  of  the  ice  cores 
from  the  Beaufort  Sea  had  a  shear  wave  velocity  that  varied  with  direction  as  shown  in  Fig. 
8.1.  The  diameter  in  the  X  direction  was  10.38  cm  and  the  shear  velocity  13  was  1730  m/s. 
The  diameter  in  the  Y  direction  was  10.44  cm  and  /3  =  1831  m/s.  Sheets  of  tiny  air  bubbles 
trapped  in  the  ice  core  caused  the  observed  variation  in  shear  wave  velocity. 

Source  and  receiver  transducers  were  shifted  vertically  by  1  cm  along  the  ice  core  and 
the  measurements  were  repeated  (Fig.  8.2).  The  diameter  in  the  X  direction  was  10.49  cm 
and  /?  =  1748  m/s.  The  diameter  in  the  Y  direction  was  10.35  cm  and  /3  =  1815  m/s. 

The  maximum  variation  in  shear  wave  velocity  along  two  orthogonal  directions  across 
the  cylindrical  core  was  about  6  %  .  while  the  diameter  changed  by  only  1.3  %  .  The 
repeatability  ot  the  results  was  better  than  5  x  lO"*.  It  is  relevant  to  mention  that  typical 
studies  of  sea  ice  shear  wave  velocity  variation  with  temperature  deal  with  velocity  fluctuations 
of  about  5  %. 

The  previous  experiments  were  repeated  on  an  outdoor-pool  sea  ice  core  sample.  The 
outdoor-pool  sample  was  produced  in  a  controlled  environment  and  yielded  more  uniform 
shear  velocity  than  the  Beaufort  sea-ice  sample.  Typical  shear  wave  velocity  variations  across 
two  orthogonal  diametrical  paths  are  shown  in  Fig,  8.3.  The  diameter  in  the  X  direction  was 
9.97  cm  and  ^  =  1749  m/s.  The  diameter  in  the  Y  direction  was  9.92  cm  and  /3  =  1771  m/s. 


8.1  References 

18.1]  J,  R.  Chamuel,  "Ultrasonic  Characterization  of  Arctic  3-D  Transient  Seismo- 
Acoustic  Phenomena,"  in  J.  Fondrk  and  K.  Dillard,  Office  of  Naval  Research 
Arctic  Sciences  Program  Research  Summaries  1992,  ONR  Report  No.  OCNR 
1125AR93-10,  May  1993,  pp.  138-140. 
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Variation  in  shear  wave  velocity  across  two  orthogonal  diametrical  paths 
measured  from  Beaufort  Sea  ice  core  sample. 
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Fig.  8.2.  Same  as  Fig. 8.1  measured  at  Z  =  Z^-\-  1  cm. 


Fig.  8.3.  Shear  wave  velocity  results  from  an  outdoor-pool  sea  ice  core  sample. 
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SECTION  9 


9  SEISMIC  PROFILING  USING  ONE  COMPRESSIONAL  TRANSDUCER  AND 
MONOSTATIC  SCATTERING  FROM  45°  OBLIQUE  CYLINDRICAL  HOLE 

9.1  Introduction 


Laboratory  ultrasonic  modeling  results  are  presented  demonstrating  the  use  of 
monostatic  scattered  waves  (2<ka<8)  from  a  45°  oblique  cylindrical  hole  in  a  solid  to 
simultaneously  determine  compressional,  shear,  and  Rayleigh  wave  velocities  (9.1).  A  vertical 
compressional  disc  transducer  placed  on  the  free  surface  of  a  solid  is  used  to  generate  and 
detect  diagonally  propagating  shear  waves  created  by  the  vertical  displacement  of  the  disc 
|9.2|.  Fig.  9.1  shows  a  diagram  by  R.  D.  Woods  |9.3|  revealing  the  diagonally  propagated 
shear  waves  generated  by  the  vertical  displacement  of  a  disc  source  on  an  elastic  half-space. 
Backscattering  results  from  45°  oblique  boundaries  in  glass,  ice  and  Plexiglas  models  are 
described.  Data  gathered  from  layered  models  show  the  effect  of  velocity  variation  with  depth 
on  ray  bending  and  monostatic  scattering.  The  proposed  oblique  hole  monostatic  reflection 
seismology  method  may  be  used  to  determine  sea  ice  shear  velocity  using  only  one 
compressional  transducer  placed  on  the  free  surface  of  the  ice.  and  a  45°  hole  bored  with  a 
coring  machine  such  as  the  PICO  45°  ice  coring  rig  currently  used  in  sea  ice  electromagnetic 
and  optical  oceanography  experiments  (W.  S.  Pegau  et  al.  |9.4|).  The  wave  propagation 
distance  between  transducer  and  retlection  point  deep  inside  the  ice  may  be  approximated  from 
the  free  surface  yielding  accurate  one-sided  velocity  measurements. 


! 


Fio.  6.19.  Distribution  of  displacement  and  energy  in  dilatational, shear. and  surface  waves 
from  a  harmonic  normal  load  on  a  half-space  for  »>  —  J.  (After  Woods  (44,  Fig.  I].) 


Fig.  9. 1  Elastic  waves  generated  by  vertical  displacement  of  disc  on  a  solid  half-space. 

(After  R.  D.  Woods,  copied  from  19. 3|).  Note  diagonally  generated  shear  wave. 
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The  basic  geometrical  configuration  of  the  proposed  method  is  diagrammed  in  Fig.  9.2. 
The  velocity  V  is  determined  from  the  relation 


where  X  is  the  horizontal  distance  between  transducer  and  hole  top,  and  t  is  the  round  trip 
propagation  time  of  waves  normally  incident  on  the  45°  oblique  hole.  The  distance  "d" 
between  the  transducer  and  the  reflection  point  along  the  oblique  discontinuity  is  equal  to 
0.707  IX.  Readers  interested  in  the  scattering  characteristics  of  cylindrical  discontinuities  may 
refer  to  papers  by  R.  M.  White  19.5]  and  Beattie  et  al.  |9.6]. 

Other  experimental  results  are  also  presented  in  this  section  on  forward  scattering  from 
a  combination  of  oblique  and  vertical  cylindrical  holes.  The  proposed  techniques  may  be 
applicable  to  NDE.  geophysics,  and  sea  ice  characterization. 


S/R 


Fig.  9.2  Geometrical  configuration  of  setup. 


9.2  Experimental  Setup  and  Results 

Ultrasonic  models  were  utilized  to  gain  insight  into  the  problem  of  monostatic 
scattering  of  compressional  and  shear  waves  inside  a  solid  by  a  45°  oblique  discontinuity  such 
as  a  cylindrical  hole.  The  first  set  of  experiments  were  conducted  on  2-D  ultrasonic  models 
as  described  below.  A  45°  wedge  was  cut  from  a  2.15  mm  thick  glass  plate.  The  glass 
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wedge  was  fixed  vertically  with  one  edge  levelled  horizontally,  and  source/receiver 
piezoelectric  compressional  transducer  was  placed  on  the  horizontal  edge  at  a  distance  X  = 
14.14  cm  from  the  45°  corner  as  shown  in  Fig.  9.3.  The  selected  distance  X  corresponded 
to  d  =  10  cm  (distance  between  the  transducer  and  the  inclined  edge  of  the  wedge).  The 
transducer  was  excited  with  a  single  broadband  pulse  and  the  received  signal  was  amplified, 
filtered,  and  displayed  on  a  digital  oscilloscope.  The  waveform  of  the  received  signal  is 
displayed  in  Fig.  9.3.  The  horizontal  time  scale  in  the  figure  is  10  us/div.  Three  reflected 
wave  components  were  detected  corresponding  to  compressional,  shear  and  Rayleigh  waves. 
The  compressional  and  shear  waves  propagated  inside  the  solid  and  reflected  from  point  A, 
while  the  Rayleigh  wave  propagated  along  the  horizontal  edge  and  reflected  from  B.  The 
velocities  of  the  plate  compressional,  shear,  and  edge  Rayleigh  waves  were  5415  m/s,  3413 
m/s,  and  3104  m/s,  respectively. 


Fig.  9.3  2-D  ultrasonic  model  simulating  reflections  of  compressional,  shear,  and 

Rayleigh  waves  from  45°  oblique  borehole  in  solid  half-space.  Transducer 
placed  along  edge  of  45°  glass  plate  wedge. 
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The  experiment  described  above  was  repeated  using  a  1  cm  thick  ice  plate.  Fig.  9.4 
shows  the  waveform  of  the  PP,  SS,  and  RR  wave  components  measured  at  X  =  7.07  cm,  9.89 
cm,  and  14.14  cm.  The  measured  velocities  were  2916  m/s,  1886  m/s,  and  1798  m/s  for  the 
plate  compressional,  shear,  and  Rayleigh  waves,  respectively.  The  horizontal  time  scale  in 
Fig.  9.4  is  20  us/div.  . 

Results  from  a  3-D  ice  model  are  presented  in  Fig.  9.5.  A  flat  hozizontal  surface  was 
prepared  on  a  large  ice  block,  a  45°  oblique  hole  was  bored,  and  the  source/receiver 
transducer  was  placed  on  the  flat  ice  surface  as  shown.  The  horizontal  distance  between  S/R 
and  the  hole  was  4.5  cm.  The  source  was  again  excited  with  a  single  broadband  pulse.  PP, 
SS,  and  RR  waves  were  observed  when  the  45°  oblique  hole  was  created  (Fig. 9. 5(a))  similar 
to  the  previous  2-D  models  results.  Tiny  air  bubbles  and  inhomogeneities  in  the  ice  block 
contributed  to  the  complexity  of  the  received  waveform.  Fig.  9.5(b)  shows  the  reference 
received  waveform  obtained  from  the  ice  block  before  introducing  the  45°  hole.  The 
horizontal  time  scale  in  Fig.  9.5  is  10  us/div.  . 


Fig.  9.4  2-D  ultrasonic  model  simulating  reflections  of  compressional,  shear,  and 

Rayleigh  waves  from  45°  oblique  borehole  in  solid  half-space.  Transducer 
placed  along  edge  of  45  °  ice  plate  wedge. 
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The  presence  of  layers  and  velocity  variation  with  depth  cause  ray  bending.  Results 
from  a  layered  2-D  ultrasonic  model  are  shown  in  Fig.  9.6  and  Fig.  9.7.  The  geometrical 
configuration  of  the  layered  model  and  the  elastic  properties  are  presented  in  Fig.  9.6(a).  The 
model  consisted  of  a  4  cm  thick  low-velocity  layer  (Cp,  =  4066  m/s,  C,,  =  2207  m/s)  over 
a  high-velocity  bottom  (Cp2  =  4857  m/s,  C,2  =  2957  m/s).  The  density  ratio  was 
approximately  equal  to  1.  The  received  waveform  is  displayed  in  Fig.  9.6(b).  The  PP  and 
SS  wave  components  arrived  at  45.3  us  and  79.28  us,  respectively,  matching  the  calculated 
results  shown  in  Fig.  9.7.  The  calculated  angle  of  incidence  at  the  layer  bottom  interface  was 
36.29°  for  the  compressional  wave  and  31.85°  for  the  shear  wave.  In  the  layered  model, 
shear  and  compressional  waves  reflected  from  different  locations  along  the  45°  oblique 
boundary. 


PP  SS  RR 

Fig.  9.5  3-D  ultrasonic  modeling  results  using  45°  oblique  hole  in  ice  block. 


Fig.  9.6  Compressional  and  shear  waves  from  2-D  layered  ultrasonic  model. 
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COMPRESSIONAL  WAVE 


SHEAR  WAVE 


45.3  MS 
PP 


T55=  79.28 


Fig.  9.7  Calculated  results  for  layered  model. 


SHEAR  WAVE  VELOCITY  (m/s)  HORIZONTAL  DISTANCE  (m) 

a  b 

Fig.  9.8  Effect  of  shear  wave  variation  with  depth  on  ray  path. 


The  effect  of  shear  wave  velocity  variation  with  depth  on  ray  bending  is  demonstrated 
in  Fig.  9.8  using  a  simplified  sea  ice  velocity  profile.  The  assumed  shear  wave  velocity 
profile  is  shown  in  Fig.  9.8(a).  The  calculated  curved  ray  path  for  the  depth-dependent  shear 
wave  velocity  model  is  plotted  in  Fig.  9.8(b).  These  results  indicate  that  shear  waves  in  the 
depth-dependent  model  would  intersect  a  45°  boundary  at  2  m  depth  instead  of  2.25  m  which 
corresponds  to  the  isotropic  model.  Further  studies  are  needed  to  characterize  layered  models. 
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Other  preliminary  experimental  results  are  also  presented  in  this  section  on  forward 
scattering  from  a  combination  of  oblique  and  vertical  cylindrical  holes  in  Plexiglas  models. 
An  interesting  configuration  is  shown  in  Fig.  9.9(a)  consisting  of  two  intersecting  45°  oblique 
cylindrical  holes  forming  a  "V"  shaped  profile  equivalent  to  a  corner  reflector.  In  an  isotropic 
model,  when  the  source  and  receiver  are  positioned  as  shown  in  Fig.  9.9(a),  the  total  travel 
time  of  the  waves  reflected  by  the  V  boundaries  is  independent  of  the  horizontal  distance  x 
between  source  and  receiver.  The  top  and  bottom  traces  in  Fig.  9.9(b)  were  recorded  with 
the  transducers  separated  by  x  =  2.7  cm  and  x  =  6.5  cm,  respectively.  Notice  the  time  of 
arrival  of  the  compressional  wave  pulse  did  not  change  as  x  was  varied. 

The  previous  configuration  shown  in  Fig.  9.9  required  two  45°  holes  and  two 
independent  source  and  receiver  transducers.  Figure  9.10  presents  a  different  configuration 
suitable  for  monostatic  applications  utilizing  one  45°  hole  intersecting  one  vertical  hole  as 
shown  and  one  source/receiver  transducer. 

9.3  Conclusions 

A  one-sided  velocity  profiling  measurement  technique  for  simultaneous  determination 
of  compressional,  shear,  and  Rayleigh  wave  velocities  using  only  one  compressional 
source/receiver  disc  transducer  and  a  45°  oblique  cylindrical  hole  was  introduced.  Ultrasonic 
modeling  results  were  presented  from  isotropic  and  layered  models  verifying  the 
implementation  of  the  proposed  method. 

Potential  advantages  of  the  proposed  seismic  profiling  technique  include:  1)  No 
instrumentation  needed  inside  the  borehole.  2)  Only  compressional  transducers  are  needed  to 
determine  both  shear  and  compressional  wave  velocities.  3)  One-sided  accurate  velocity 
measurements  may  be  achieved.  4)  Dual  use  of  hardware  and  45°oblique  holes  utilized  in 
sea  ice  optical  oceanography  experiments  |9.4J. 


b 


a 


Fig.  9.9  V  configuration  in  Plexiglas  model  with  45°  drilled  holes.  Notice  arrival  time 
of  compressional  pulse  is  independent  of  x. 
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AIR  FILLED  BOREHOLE 


WATER  FILLED 


Fig.  9.10  Modified  V  configuration  for  monostatic  applications.  Filling  the  6.3  mm  dia 
cylindrical  cavities  with  water  drastically  reduced  the  reflected  compressional 
wave. 
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SECTION  10 


10  SEA  ICE  FINGER-RAFTING  AND  ICE  THICKNESS 
10.1  Introduction 

Sea  ice  finger-rafting  is  a  periodic  type  of  ice  structure  in  the  form  of  periodic  parallel 
fractures  often  perpendicular  to  the  edges  of  two  adjacent  floating  thin  ice  plates  with 
alternating  rectilinear  overthrusts  and  underthrusts  (lO.lJ.  The  fracture  pattern  is  at  times  so 
regular  that  Weeks  and  Anderson  [10. 1]  mentioned  that  many  of  the  rectilinear  thrusts  appear 
to  have  been  laid  out  by  a  draftsman.  Fig.  10. 1  compares  the  finger-rafting  diagram  of  Weeks 
and  Anderson  with  that  of  Homes  and  Worthington  [10.21  as  reported  by  Dunbar  [10.3]. 
Weeks  and  Anderson  described  the  formation  process  as  an  advancing  ice  sheet  interacting 
with  a  stationary  fast  ice  edge,  where  fracturing  occurs  only  in  the  moving  ice  sheet.  Homes 
and  Worthington  described  finger-rafting  as  fracture  and  interlocking  of  two  ice  edges  in 
refrozen  leads.  Weber  [10.4|  and  Green  [10.5|  observed  finger-rafting  in  lakes. 

Weeks  and  Anderson  observed  that  there  is  no  apparent  correlation  between  the  ice 
plate  thickness  and  the  width  or  length  of  the  finger-rafting  features.  Although  Weeks  and 
Anderson  did  not  observe  the  actual  process  of  finger-rafting  formation,  they  provided  an 
explanation  suggesting  that  the  extremely  flexible  newly  formed  sea  ice  transmits  water  waves. 
The  ragged  edge  of  an  advancing  ice  sheet  is  deposited  on  top  of  another  ice  sheet  by  the  crest 
of  a  water  wave.  The  deposited  ice  sheet  tears,  and  part  of  the  collapsed  portion  slides  under 
the  fast  ice  sheet  during  the  following  wave  trough.  The  process  advances  to  the  adjacent  ice 
edges  segments  under  the  continued  action  of  wind  and  water  wave  pushing  the  offshore  sheet 
into  the  shore  fast  ice  [10.11. 

Dunbar  [10.31  suggested  an  explanation  for  the  formation  of  finger  rafting  when  the 
water  waves  move  at  an  angle  with  respect  to  the  wind  direction  and  are  obliquely  incident  on 
the  adjacent  parallel  edges  of  the  two  ice  plates.  Dunbar  mentioned  that  the  finger-rafting 
width  depends  on  the  wavelength  and  angle  of  incidence  and  will  bear  no  relation  to  the  ice 


Wttks  and  Anderson 


Holtnes  and  Worlhir^gton 


Fig.  10.1  Finger-rafting  diagram  by  Weeks  and  Anderson  [10.11  compared  with  a 
diagram  by  Homes  and  Worthington  [10.21  as  reported  by  Dunbar  110.31. 
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FREQUENCY  (Hz) 

Fig.  10.2  Phase  velocities  of  water  gravity  wave  and  flexural  wave  in  floating  sea  ice 
plate  (5  cm  thick)  plotted  as  function  of  frequency. 

thickness.  One  puzzling  detail  in  Dunbar’s  paper  is  in  his  Figure  10,  where  the  fractures  were 
sketched  perpendicular  to  the  plate  edges  independent  of  the  water  wave  angle  of  incidence. 
Dunbar’s  explanation  of  the  formation  process  depends  on  a  propagating  process  along  the  ice 
plate  edge  originating  from  one  corner  of  the  floating  plate  and  progressing  along  the  edge 
following  the  water  wave.  Physically  speaking,  the  edge  of  the  plate  cannot  overlap  the 
second  ice  plate  unless  the  whole  plate  moves  simultaneously  along  each  point  along  the  edge. 
The  ice  plate  is  flexible  and  can  bend  vertically,  however,  it  may  not  be  possible  to  lift  one 
corner  of  the  plate  and  at  the  same  time  elongate  the  plate  to  have  it  deposited  on  top  of  the 
second  plate. 

10.2  Results 

The  only  known  quantitative  data  on  finger-rafting  is  a  single  point  obtained  by  J. 
Morrison  1 10.7)  observed  in  5  cm  thick  sea  ice  where  the  measured  periodic  finger  width  was 
5  m.  R.  S.  Pritchard  [10.7]  presented  a  theoretical  relationship  between  finger  width  and  ice 
plate  thickness.  Pritchard’s  calculated  finger-rafting-width  results  were  smaller  than 
Morrison’s  single  data  point  by  almost  a  factor  of  5.  Other  qualitative  Arctic  field 
observations  by  T.  B.  Curtin  (1990-1993)  indicated  that  finger-rafting  occurs  over  a  vary  wide 
range  of  ice  thicknesses.  For  an  ice  thickness  of  about  1  m,  Curtin  observed  finger-rafting- 
width  between  50-100  m. 

New  theoretical  results  were  obtained  by  the  author  [10.8]  relating  finger-rafting-width 
to  ice  plate  thickness.  The  results  were  based  on  the  assumption  that  when  the  flexural  wave 
speed  in  a  floating  plate  (neglecting  gravity)  is  equal  to  the  speed  of  a  water  gravity  wave 
component,  maximum  coupling  between  the  gravity  wave  and  flexural  wave  occurs.  The 
theoretical  results  were  obtained  by  calculating  the  intersection  point  between  the  flexural  wave 
phase  velocity  dispersion  curve  of  a  floating  ice  plate  and  the  phase  velocity  dispersion  curve 
of  water  gravity  waves  calculated  at  different  thicknesses  (Fig.  10.2).  The  gravity  wave  phase 
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velocity  plot  was  obtained  from  the  expression; 


v(0  =  g/27rf 

where  v  is  the  phase  velocity,  g  is  gravity,  and  f  is  frequency. 

The  flexural  wave  phase  velocity  was  calculated  from  the  characteristic  equation  [10.9]; 


pv\  1-. 


GiV)  *  D(V)  2p\  1-- 


where  c  =  water  compressional  wave  speed 
a  =  ice  compressional  wave  speed 
/3  =  ice  shear  wave  speed 
p  =  ratio  of  water  density/ice  density 


G{V)  = 


coth  I 


■m  (i-z!) 

^  V 


2 


coth 


(1--^)) 


D(V)  = 


tanh  ( 


t/2  2 


(1-^)  tanh(il^,  (1-i?)) 


The  intersection  point  between  the  flexural  wave  phase  dispersion  curve  and  the  gravity 
wave  dispersion  curve  yielded  the  desired  phase  velocity  and  frequency  from  which  the 
wavelength  was  calculated.  Finger-rafting  occurs  at  half  this  wavelength.  The  calculated 
finger-rafting  width  was  plotted  as  function  of  ice  plate  thickness  in  Fig.  10.3(a).  An 
expanded  plot  is  presented  in  Fig.  10.3(b)  indicating  that  for  5  cm  thick  ice,  the  calculated 
finger-width  is  5.3  m.  Morison’s  data  point  and  Pritchard’s  theoretical  results  were  also 
plotted  in  Fig.  10.3(b)  for  comparison.  The  parameters  used  in  the  calculations  were  as 
follows;  ice  compressional  wave  velocity  =3130  m/s,  ice  shear  wave  velocity  =  1670  m/s, 
water/ice  density  ratio  =  1.025,  water  compressional  wave  velocity  =  1440  m/s,  gravity  = 
9.80665  m/s\  Numerical  results  indicated  that  increasing  the  elastic  parameters  of  the  ice  to 
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FINGER  RAFTING  WIDTH  (H) 


ICE  PLATE  THICKNESS  (M) 


ICE  PLATE  THICKNESS  (M) 


Fig.  10.3  A)  Calculated  fmger-rafting-width  plotted  as  function  of  ice  thickness.  Solid 
line  with  open  squares  (Chamuel)  and  dashed  line  (Pritchard).  B)  Expanded 
scale  showing  LEADEX  Arctic  field  data  point  by  Morrison  (solid  black  circle). 
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match  hard-ice  had  little  effect  on  the  predicted  finger-width  values.  For  example,  the 
predicted  finger-rafting  width  for  1  m  thick  ice  changed  only  from  50  m  to  51  m  as  the  above 
ice  properties  were  replaced  with  corresponding  hard-ice  elastic  properties. 

An  empirical  equation  was  obtained  to  approximate  the  theoretical  results  relating  the 
calculated  finger-width  W  to  ice  thickness  H; 

W  =  98  H  ° 

Fig.  10.4  compares  the  empirical  and  theoretical  results.  Finger-rafting  is  difficult  to  simulate 
in  the  laboratory  using  a  scale  model.  Ice  fracture  in  the  Arctic  depends  on  the  unique 
properties  of  sea  ice  at  the  natural  scale. 


10.3  Conclusions 

A  preliminary  study  was  conducted  to  provide  physical  insight  into  the  formation 
process  relating  finger-rafting  width  to  ice  plate  thickness.  Theoretical  results  were  obtained 
based  on  the  assumption  that  maximum  coupling  between  flexural  waves  and  water  gravity 
waves  occur  when  their  phase  speeds  are  identical.  The  theoretical  results  were  compared  to 
other  theoretical  and  field  results. 

Fundamental  theoretical,  field,  and  laboratory  studies  on  sea  ice  finger-rafting  are 
needed  to  complement  the  preliminary  investigations  conducted  so  far.  Future  laboratory 
experimental  studies  will  characterize  the  action  of  coupled  flexural  waves  and  gravity  waves 
in  adjacent  finite  flexible  floating  plates. 


ICE  PLATE  THICKNESS  (M) 


Fig.  10.4  Comparison  of  empirical  (solid  line)  and  theoretical  (open  squares)  results. 
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SUMMARY 


This  report  is  the  third  in  a  series  of  publications  by  the  author  [5.5,  5.13],  presenting 
a  new  collection  of  ultrasonic  modeling  findings  obtained  under  controlled  conditions  on 
various  fundamental  liquid/solid  seismoacoustic  wave  phenomena.  Quantitative  ultrasonic 
modeling  results  were  described  contributing  to  the  physical  understanding  of  different 
fundamental  problems  related  to  low-frequency  seismoacoustic  wave  phenomena  resulting  from 
the  interaction  of  underwater  acoustic  waves  with  penetrable  heterogeneous  solid  boundaries. 
The  new  findings  were  summarized  in  ten  independent  sections  covering  a  wide  range  of 
problems: 


1)  Backscattering  of  Scholte  waves  and  grazing  underwater  acoustic  waves  by 
a  trench  at  a  liquid/solid  interface. 

2)  Explanation  for  observed  underwater  acoustic  backscatter  from  the  Canada 
Basin  near  73.2°N,  I39°W  where  recent  bathymetric  charts  do  not  reveal  the 
existence  of  large  bathymetric  highs  nearby. 

3)  Effect  of  a  thin  low-velocity  layer  underneath  a  floating  plate  on  flexural 
wave  dispersion. 

4)  Effect  of  viscous  waves  from  suspended  particles  in  a  thin  layer  under  sea 
ice  on  low-frequency  underwater  acoustic  wave  reflectivity. 

5)  Comparison  of  experimental  and  theoretical  results  on  the  dispersion  of 
antisymmetric  edge  waves  along  apex  of  truncated  elastic  wedge. 

6)  Wave  propagation  along  apex  of  free  and  immersed  elastic  wedges  with 
range-dependent  apex  angle  and  cross  section. 

7)  Shear  wave  coupling  to  ice  cores. 

8)  Ultrasonic  characterization  of  sea  ice  cores. 

9)  Monostatic  seismic  profiling  using  scattering  from  45°  oblique  cylindrical 
hole  in  ice  plate  using  one  compressional  transducer. 

10)  Sea  ice  finger-rafting  and  ice  thickness. 

The  generic  findings  presented  are  applicable  to  the  fields  of  underwater  acoustics, 
acoustical  oceanography,  seismology,  nondestructive  testing,  acoustoelecronic  sensors,  acoustic 
imaging,  and  noise  reduction. 
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